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Abstract 

This paper is organized in two parts. The first part (Sections 2-5) is dedicated to the 
theory at T = and contains a pedagogical review of some fundamental aspects related 
with the chiral symmetries of QCD, the U(l) problem and its solution proposed by 'tHooft, 
Witten and Veneziano. In the second part (Sections 6-14) we discuss the role of the U(l) 
axial symmetry for the phase structure of QCD at finite temperature. One expects that, 
above a certain critical temperature, also the U(l) axial symmetry will be restored. We 
will try to see if this transition has (or has not) anything to do with the usual chiral 
transition: various possible scenarios are discussed. In particular, we analyse a scenario 
in which the U(l) axial symmetry is still broken above the chiral transition. We will show 
that this scenario can be consistently reproduced in the full respect of the relevant QCD 
Ward Identities and also using an effective Lagrangian model. A new order parameter is 
introduced for the U(l) axial symmetry. 



1. Introduction 



It is generally believed that a phase transition which occurs in QCD at a finite temper- 
ature is the restoration of the spontaneously broken SU(L) <S> SU(L) chiral symmetry in 
association with L massless quarks. At zero temperature the chiral symmetry is broken 
spontaneously by the condensation of qq pairs and the L 2 — 1 J p = mesons are just 
the Nambu-Goldstone (NG) bosons associated with this breaking fl], ||, ||, |], |]. At high 
temperatures the thermal energy breaks up the qq condensate, leading to the restoration 
of chiral symmetry. We expect that this property not only holds for massless quarks but 
also continues for a small mass region. The order parameter for the chiral symmetry 
breaking is apparently (qq) = J2i=i(QiQi) '■ the chiral symmetry breaking corresponds to 
the non-vanishing of (qq) in the chiral limit sup(mj) — > 0. From lattice determinations of 
the chiral order parameter (qq) one knows that the SU(L)<g>SU(L) chiral phase transition 
temperature T ch , defined as the temperature at which the chiral condensate (qq) goes to 
zero (in the chiral limit sup(mj) — > 0), is nearly equal to the deconfining temperature T c 
(see, e.g., Ref. ||). But this is not the whole story: QCD possesses not only an approxi- 
mate SU(L) ® SU(L) chiral symmetry, for L light quark flavours, but also a U(l) axial 
symmetry (at least at the classical level) [0, §]. The role of the U(l) symmetry for the 
finite temperature phase structure has been so far not well studied and it is still an open 
question of hadronic physics whether the fate of the U(l) chiral symmetry of QCD has 
or has not something to do with the fate of the SU(L) ® SU(L) chiral symmetry. In the 
following Sections we will try to answer these questions: 

• At which temperature is the U(l) axial symmetry restored? (if such a critical 
temperature does exist!) 

• Does this temperature coincide with the deconfinement temperature and with the 
temperature at which the SU(L) chiral symmetry is restored? 

In the "Witten-Veneziano mechanism" ||, [K| for the resolution of the U(l) problem, a 
fundamental role is played by the so-called "topological susceptibility" in a QCD without 
quarks, i.e., in a pure Yang-Mills (YM) theory, in the large-iV c limit (N c being the number 
of colours): 




(1.1) 



2 



2 

where Q(x) = ^2^ upa F^ v F^ a is the so-called "topological charge density" . This quantity 
enters into the expression for the squared mass of the rf: nil, = ^kr, where L is the number 
of light quark flavours taken into account in the chiral limit. Therefore, in order to study 
the role of the U(l) axial symmetry for the full theory at non-zero temperatures, one 
should consider the YM topological susceptibility A(T) at a given temperature T, formally 
defined as a large-iV c limit of a certain expectation value (■ ■ -)t i n the full theory at the 
temperature T fllTf : 

A(T) = lim lim l-i [ d 4 xe ikx (TQ(x)Q(0)) T ) . (1.2) 

fc^O N c —*oo I J J 

In other words, ACT) = lim A (k,T) = A (0,T), where the quantity 

k— >0 

A (k,T) = lim l-i [ d 4 xe ikx (TQ(x)Q(0)} T \ , (1.3) 

at four-momentum k and physical temperature T, is nothing but the leading-order term 
in the 1/N C expansion of the corresponding quantity 

X (k, T) = -il d 4 xe* kx (TQ(x)Q(0)) T , (1.4) 



of the full theory, at the same four-momentum k and physical temperature T. [From 
the general theory of 1/N C expansion it is known that Ao(k,T) is of order 0(1), while 
successive terms coming from including fermions, are of order 0(1 /N c ).] It is in this 
sense that one can talk about the behaviour of A(T) above or below the chiral transition 
temperature T c h. 

The problem of studying the behaviour of A(T) as a function of the temperature T 
was first addressed, in lattice QCD, in Refs. [0, [13], [L4]]. Recent lattice results [15[ 



(obtained for the 5*^7(3) pure-gauge theory) show that the YM topological susceptibility 
A(T) is approximately constant up to the critical temperature T c ~ T c h, it has a sharp 
decrease above the transition, but it remains different from zero up to ~ 1.2 T c . In the 
Witten-Veneziano mechanism |], OT, a (no matter how small!) value different from zero 
for A is related to the breaking of the U(l) axial symmetry, since it implies the existence 
of a pseudo-Goldstone particle with the same quantum numbers of the rf (see also Ref. 
16|). Therefore, the available lattice results for the topological susceptibility show that 



the U(l) chiral symmetry is restored at a temperature Tu(i) greater than T c h- 
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Another way to address the same question is to look at the behaviour at non-zero 
temperatures of the susceptibilities related to the propagators for the following meson 
channels [I7| (we consider for simplicity the case of L = 2 light flavours): the isoscalar 



1 = scalar channel a (also known as fo in the modern language of hadron spectroscopy), 
interpolated by the operator O a = qq; the isovector 1=1 scalar channel 5 (also known 
as do), interpolated by the operator 0$ = q\q\ the isovector 1=1 pseudoscalar channel 
7T, interpolated by the operator O n = i(pi$\q\ the isoscalar 1 = pseudoscalar channel 
T)', interpolated by the operator O v ' = iq^q. Under SU(2)a transformations, a is mixed 
with 7r: thus the restoration of this symmetry at T c h requires identical correlators for 
these two channels. Another SU(2) chiral multiplet is (8,rf). On the contrary, under 
the U(1)a transformations, 7r is mixed with 5: so, a "practical restoration" of the U(l) 
axial symmetry should imply that these two channels become degenerate, with identical 
correlators. Another U(l) chiral multiplet is (a, rj'). (Clearly, if both chiral symmetries 
are restored, then all tt, rj', a and 5 correlators should become the same.) In practice, one 
can construct, for each meson channel /, the corresponding chiral susceptibility 

Xf = Jd 4 x(O f (x)O}(0)), (1.5) 

and then define two order parameters: Xsu(2)®su(2) = Xa - Xtt, and xu{i) = Xs - X-n- 
If an order parameter is non-zero in the chiral limit, then the corresponding symmetry 
is broken. Present lattice data for these quantities seem to indicate that the U(l) axial 
symmetry is still broken above T C h, up to ~ 1.2 T c h, where the 5-tc splitting is small but 



still different from zero |T^, [19|, |2(| . In terms of the left-handed and right-handed quark 
fields [qL,R = |(1 ± 7s)<?], one has the following expression for the difference between the 
correlators for the S + and ir + channels: 

V u(1) (x) = (O 5+ (x)Oj + (0)) - (O^(x)OU(0)) 

= 2 \{u R d L {x) ■ d R u L {0)) + (u L d R (x) ■ d L u R {0))] . (1.6) 

(The integral of this quantity, / d 4 xVu^(x), is just equal to the U(l) chiral parameter 
Xu(i) — Xs ~ X-k-) What happens below and above T ch l Below T ch , in the chiral limit 
sup(mj) — > 0, the left-handed and right-handed components of a given light quark flavour 
(up or down, in our case with L = 2) can be connected through the quark condensate, 
giving rise to a non-zero contribution to the quantity T>u(i)(x) in Eq. (|1.6| ) (i.e., to the 
quantity Xu{i))- But above T c h the quark condensate is zero: so, how can the quantity 
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T>xj{\){x) (i.e., the quantity Xu(i)) be different from zero also above T c h, as indicated by 
present lattice data? The only possibility in order to solve this puzzle seems to be that 
of requiring the existence of a genuine four-fermion local condensate, which is an order 
parameter for the U(l) axial symmetry and which remains different from zero also above 
T ch . This will be discussed in Section 7. 

The paper is organized in two parts. The first part (Sections 2-5) is dedicated to 
the theory at T = and contains a pedagogical review of some fundamental aspects 
related with the chiral symmetries of QCD (Section 2), the U(l) problem and its solution 
proposed by 'tHooft (Section 3), the Witten's mechanism (Section 4) and the Veneziano's 
mechanism (Section 5). In the second part (Sections 6-14) we discuss the role of the U(l) 
axial symmetry for the phase structure of QCD at finite temperature. One expects that, 
above a certain critical temperature also the U(l) axial symmetry will be restored. 

We will try to see if this transition has (or has not) anything to do with the usual chiral 
transition: various possible scenarios are discussed in Section 6. In particular, we analyse 
a scenario in which the U(l) axial symmetry is still broken above the chiral transition 



2T|, |22|, 23"|. In Section 7 a new order parameter is introduced for the U(l) axial symmetry. 



In Section 8 we make an analysis of the relevant QCD Ward Identities (WI) in the range of 
temperatures between T C h and Tjj^ . The saturation of the relevant QCD Ward Identities 
is obtained in a 1/JV C expansion. A new 2L-fermion effective vertex is introduced to take 
care of the new order parameter for the U(l) chiral symmetry alone. In Sections 8 and 
9 we introduce a new meson field associated with the new U(l) chiral order parameter 
and then, in Section 9, we study the form of a new effective Lagrangian including, in 
addition to the usual quark-ant i- quark meson fields Uij, also the new 2L-fermion field X. 
In Sections 10 and 11 we study the mass spectrum of the theory respectively above and 
below the SU (L) <g> SU (L) chiral phase transition. In Section 12 we apply the results of 
Section 11 to the "real-world' 1 case (at T < T c h), where there are L = 3 light flavours, 
named u, d and s: we derive a generalized Witten-Veneziano formula for the rj' mass. 
In Section 13 we show that the above results are "model-independent", i.e., independent 
on additional terms in the Lagrangian. In Section 14 we analyse the consequences of our 
theoretical model on the chiral condensate (qq) and on the topological susceptibility x, 
in the full theory with quarks. In particular we are interested in their dependence on the 
light quark masses (something about this is also said in the Appendix). In Section 14 we 
also discuss the ^-dependence for some relevant quantities that we have found. Finally, 
the conclusions and an outlook are given in Section 15. 
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2. The chiral symmetries of QCD 



The parameter Aqcd ~ 0.5 GeV (in the MS renormalization scheme), associated with 
asymptotic freedom, defines a high-momentum (or short-distance) regime [k 2 / K 2 q CD ^> 
1) in which quarks and gluons can be treated as weakly interacting particles in perturba- 
tive QCD. At the other end of the scale (k 2 / Aq CD < 1) is low-energy hadronic physics, 
in which the interacting units are not individual quarks and gluons, but hadrons. So far 
we cannot solve QCD in this domain. However, from the vast amount of experimental 
information available, some highly fruitful ideas had been developed long before QCD 
was invented. Translated into expected properties of QCD, some of these ideas become 
concrete statements that are much easier to comprehend. 

Perhaps the most relevant example of this is represented by the chiral symmetries of 
strong interactions. They were known since the 60's, when the only theoretical instru- 
ments available for their analysis were current algebra and a rudimental quark model 
@. Translated into the QCD language, they become a much more clear and well-defined 
subject. Let us consider, in fact, the QCD Lagrangian: 

C-QCD = -\F; v F a ^ + qi^D.q - qMq, (2.1) 
where M is the quark mass matrix: 



M 



m d 

m s 



[2.2) 



Since the quark masses are different from one another, the symmetry group of this theory 
is G = U(l) u (g) C7 <8> • • ■ There is one conserved charge for every one of the quark 
flavours, the corresponding conserved currents being u'j^u, dj^al, etc. If we take L of 
these flavours to have zero masses (the physical meaning of this will be explained below 
0. Q|) the symmetry group of the theory becomes larger. The Lagrangian is invariant 
under independent rotations of the left- and right-handed quark fields: 

Ql -> V L q L , q R -> V R q R , (2.3) 
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where q stands for the vector (gi, q 2 , . . . , qi), while qi = |(1 + 75)9 and q R = |(1 — 75)?, 
with 75 = — i7°7 1 7 2 7 3 , denote respectively the left-handed and the right-handed quark 
fields; Vl and Vr are arbitrary L x L unitary matrices. So the symmetry group is: 

G = SU(L) L ® SU(L) R ® C/(1) L ® (2.4) 

We do not consider any more the part of the Lagrangian depending on the "heavtf quark 
fields. 

There are two ways in which a continuous global symmetry of the Lagrangian may be 
realized in the vacuum: 

• The Wigner-Weyl realization, in which all the generators of the symmetry group 
annihilate the vacuum. Then the physical states can be classified according to the 
irreducible unitary representations of the group. 

• The Nambu-Goldstone realization, in which some generators do not annihilate the 
vacuum. Then Goldstone theorem |5j states that, for each generator which does 
not annihilate the vacuum, there exists a spin-zero massless particle (Goldstone 
boson), with the quantum numbers of this generator and transforming according to 
the irreducible representations of the unbroken subgroup. 

In nature U(1)l+r is realized in the Wigner-Weyl way and manifests itself directly as 
baryon-number conservation. On the contrary SU(L)l <8> SU(L)r cannot be realized in 
such a way, because this would imply that each hadron multiplet is accompanied by a 
mirror multiplet of the same mass, but with opposite parity. No hint of this can be found 
in the hadronic spectrum. For example, there is not even an approximate mirror image of 
the nucleon iso-doublet. Thus one generally assumes that this symmetry is spontaneously 
broken down to: 

H = SU(L) L+R , (2.5) 

the vacuum being invariant only under the subgroup generated by the vector currents 
g7 M |A a g. As a consequence, the spectrum of QCD (with L massless quarks) must contain 
L 2 — 1 Goldstone bosons. Their quantum numbers can be read off from those of the states 
Qa |0), which result if the axial charge operators: 

Qi = J d 3 *A° a (x) ; A£ = qY^-Kq, (2.6) 
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are applied to the vacuum: J p = 0~. In reality, of course, the symmetry is explicitly 
broken by the quark-mass term in the Lagrangian, so that we do not observe in the 
hadronic spectrum the presence of massless Goldstone bosons. Nevertheless, the eight 
lightest hadrons (7r =t , 7r°, K , K°, K°, if) are pseudoscalars and three of them (tt , 7r°) are 
particularly light. This pattern is to be expected if m u , m^, m s happen to be small when 
compared to the typical mass-scale of strong interactions, i.e, Aq C d ~ 0.5 GeV (in the 
MS renormalization scheme). In such a way the Lagrangian is approximately invariant 
under SU(3)l ® SU(3)r and the spontaneous breakdown of this approximate symmetry 
generates 3 2 — 1 = 8 approximate massless particles. Moreover m u and must be small 
compared to m s , so that the group SU{2)l ® SU{2)r is a much better symmetry of the 
Lagrangian. This is the reason why 2 2 — 1 = 3 of the pseudoscalar mesons are particularly 

light si. 

In other words, the piece of the Lagrangian containing the lightest quark masses: 

5Cqcq = —{m u uu + m d dd + m s ss), (2.7) 



acts as a perturbative term ( "chiral perturbation theory") with respect to the remnant 
piece C^qcd °f the Lagrangian: Cqcd * s invariant under SU(3)l ® SU(3)r. On the 
contrary, the mass parameters for the other quarks c, b, t should all be much larger than 
m u , m d , m s (m c , m b , m t > Aq C d), because we see no trace of flavour S'f/(4) or higher 
symmetries in the hadronic spectrum. 

What is the dynamical reason for which SU(3)l ® SU(3)r is broken in the Goldstone 
mode? We do not have an answer. In the non-renormalizable pre-QCD model of Nambu 
and Jona-Lasinio [23 1, the cause of spontaneous symmetry breakdown is a direct nucleon- 
nucleon attraction built into the model, in analogy with the effective electron-electron 
attraction responsible for the formation of Cooper pairs in the theory of superconductivity. 
(In this theory the Goldstone boson is "eated" by the electromagnetic field through the 
Higgs mechanism and the photon becomes massive in a superconductor: this is called the 
"Meissner effect".) A direct quark-quark interaction is ruled out by renormalizability in 
QCD; but an effective interaction could arise just as the effective electron-electron attrac- 
tion in superconductivity arises from the more fundamental electron-phonon interaction. 
Maybe instantons play such a role in an effective quark-quark interaction f25fl. 

So far we have discussed what happens in nature of the symmetry SU{3)l® SU{3)r® 
U(1)l+r of the QCD Lagrangian. However, the symmetry group (|2.4|) also contains the 
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group U(1)l-r, which is usually called the "U(l) chiral group" or U(1)a (A stands for 
axial; the U(1)a group transformation is: q — *> e~ iar<5 q, q — > qe~ tai5 ). The question of how 
U{1)a is manifested presents a puzzle, which is known as the "[1(1) problem". 

If the symmetry U(1)a were manifested directly in the Wigner-Weyl mode, then in 
the chiral limit all massless hadrons would have a massless partner of opposite parity: 
but this is not the case in nature. On the other hand, if it were broken in the Goldstone 
mode, then there should be an I = pseudoscalar Goldstone boson, whose perturbed 
state should have about the same mass as the pion. Using chiral perturbation theory, 
Weinberg [7[] estimated its mass to be less than y/3m n . Among the known hadrons, the 
only candidates with the right quantum numbers are 7/(549) and 7/(958). Both violate 
the Weinberg bound: besides 7/(549) has already been claimed by the pion octet. 

The U(l) puzzle is: where is the extra Goldstone boson? 



3. The solution of the U(l) problem 

'tHooft || removed the puzzle by observing that the U(1)a is not a symmetry of the 
theory at the quantum level, for the presence of instanton effects. We will not go into the 
detailed analysis but merely mention some relevant points. 

The £7(1) axial current = J2i=i Qil^lh^i (with 75 = — 77°7 1 7 2 7 3 ) is not conserved, 
due to a QCD axial anomaly (the so-called "Adler-Bell-Jackiw (ABJ) anomaly" j26|): 

= 2LQ(x), (3.1) 

where Q(x) = -^2^ lJpa F^F^ (with e 0123 = 1) is the so-called "topological charge den- 
sity", and L is the number of light quark flavours taken into account in the chiral limit 
(L = 2 or L = 3). It is known that Q(x) is a total 4-divergence: 

Q(x) = PK ll (x) ; K, = J^s^AZ (d^Al + ~^ c A^f) ■ (3-2) 

(fabc are the group constants of SU(3) co i our ). The integral of Q(x) over all space-time 
is called the "topological charge". We can define a conserved but non-gauge-invariant 
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current: 



j(L) 



J, 



5,/j 



2LK„. 



(3.3) 



The generator of the U(1)a symmetry may be taken to be: 



Qb — I ^ 3x ^5.0 



L 

E 

i=l 



(3.4) 



In the Abelian case is g au g e invariant, because of the absence of topological charge. 
This is no longer true for the non-Abelian case, and Q5 is not a physical quantity. 



Moreover, Q\' = J d 3 x J5 is not conserved, because of the existence of instantons [57[ . 
To see this, one can integrate Eq. ( |3.1| ) over the Euclidean four-space. The result can be 
presented in the form: 
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7f 



2Lg[F], 



where: 



q[F] 



9' 



647T 5 



r/ 4 rr F a F a 



(3.5) 



(3.6) 



is the topological charge, a functional of the gauge field For F" corresponding to 
one instanton, g[F] = 1 (see Ref. [|]). Therefore, in this case, the boundary values of 
in Euclidean time differ by: 



AQi L) 



2Lq[F). 



(3.7) 



This can be attributed to the fact that an instanton interpolates (in Euclidean time) 
between two gauge-field configurations differing by one unit of topological charge. Thus 
there seems to be no reason to expect U(1)a to be a symmetry. 

However, the work of 'tHooft did remove the puzzle in its original form, but generated 
new questions. In particular: what becomes of the would-be Goldstone boson? Or, in 
other terms: why is the mass of the ?/(958) (the would-be Goldstone boson of the U{1)a) 
so large compared with the masses of the octet mesons? Witten |J and Veneziano [T(J 
proposed in 1979 a mechanism to answer these two questions and to give a reliable estimate 
of the f]' mass. This is now universally known as the "Witten-Veneziano mechanism" . 
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4. Witten's mechanism 



Let us briefly review which was the idea of E. Witten || to explain the large mass of the 
f]'. Starting from the resolution of the U(l) problem given by 'tHooft in 1976 || (which 
we have discussed in Section 3), one may try to explain the large mass of the rj' by saying 
that the Goldstone boson receives a mass as a result of the anomaly. This statement has 
not a clear meaning, however, because there is apparently no way to turn off the anomaly, 
and therefore no sense in talking of the Goldstone boson that would have existed if there 
had been no anomaly. Witten observed, however, that in a certain sense it is possible to 
vary a parameter in QCD so as to turn off the anomaly: this parameter is the number N c 
of colours of the theory. 

The U(l) problem is then reconsidered from the point of view of the 1/N C expansion 
(for a detailed analysis of this expansion see Refs. |28| and ||29|). It is for large N c that 



the anomaly turns off and it is in the sense of the 1/N C expansion that one can say that 
the anomaly gives a mass to a U(l) particle that otherwise would have been massless. 

There is also a more practical way to discuss the £7(1) problem. Isgur |JD| and De 
Rujula, Georgi and Glashow [|3TJj suggested that in the rj', which is an £77(3) singlet, the 
quark and antiquark can annihilate into gluons (this annihilation channel is absent for 
the non-singlet mesons). These quark-antiquark annihilation diagrams split the rf from 
the light pseudoscalars. Of course this is not true at a given finite order of perturbation 
theory, because in finite orders of perturbation theory we have neither an rj' nor a 71". Yet 
in the context of the 1/N C expansion, one can easily make the summation of an infinite 
class of Feynman diagrams, so as to generate bound states [3"2]] . In the large-iV c limit 
the rf and it are degenerate; this degeneracy is lifted by quark-antiquark annihilation 
diagrams that are suppressed by one power of 1/N C . 

To see how the rf gets a mass from the anomaly, we have to consider the Fourier 
transformed of the correlation function of the topological charge Q(x): 

X (k) = -ij d i xe ikx (TQ(x)Q(0)}. (4.1) 

The x(k) at zero momentum is called the topological susceptibility, which in the literature 
is usually indicated with A 01 x : 

x = -if d 4 x(TQ(x)Q(0)). (4.2) 
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It is well known, from a direct analysis of the QCD Ward Identities [33], that x is zero 
in the theory with massless quarks. This can also be seen if we write x hi the form (here 
VT = J d A x is an infinite four- volume which must be factorized): 



1 i d 2 Z g [r],T]} _ 
X VTZ e [v,fj] d9 2 1 } 

where Z$[r], fj] is the partition function of the theory including the so-called #-term (which 
in quarkless QCD labels different ^-worlds that are physically distinct): 



z e[v,v]= J [dA f 
x exp (i J d 4 x 



- l -F; u F^ + g(i 7 "£>„ - M)q + f}q + <lV + 0Q{x) 



(4.4) 



If the theory contains massless quarks [i.e., M = in Eq. ( f4.4|) l, then the partition 
function comes out to be independent of 9, being (see also Ref. [|34]|): 

zf\r 1 ^] = zfl 2La [ri\r 1 % (4-5) 

where rj = e~ iai5 r] and fj' = fje~ lOLTo (with 75 = — i7°7 1 7 2 7 3 ), and L is the number 
of massless flavours. Thus one can always reduce 9 to zero by making a global chiral 
transformation on the fermion sources, with a = —8/2L. Such a chiral transformation 
does not change the physics, because all Green's function are evaluated in the sourceless 
limit. Thus all worlds are physically equivalent to one another, when the theory contains 
massless quarks. In particular, the derivative in ( |4.4| ) vanishes and so x is equal to zero. 

Let us assume that, to leading order in 1/N C , x{k) is non-vanishing at k — 0, in the 
world without quarks. When massless quarks are included, x{k) must vanish at k = 0; 
yet the massless quarks are effects of order 1/N C , since, as we know from 'tHooft's work 
p8| , each quark loop is suppressed by a factor of 1/N C . x(h) nas the form: 



X (k) = A (k) + A x {k) + A 2 {k) + . . . , (4.6) 

where A Q (k) is the sum of all diagrams without quark loops, A\{k) is the sum of all 
diagrams with one quark loop, A 2 (k) is the sum of all diagrams with two quark loops, 
and so on. If one computes Ao(k) in perturbation theory, one finds that it is of order 
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g A N 2 = 1, since it receives contributions from each of the N 2 — 1 gluon degrees of freedom 
moreover each Q(x) takes a factor g 2 which is of order 1/N C , since, as 'tHooft showed |28 



to have a smooth large-iV c limit one should define the coupling constant to be g = g / \HST C , 
where go is kept fixed as N c becomes large. Yet Ai(k) is of order 1/N C , A 2 (k) is of order 
l/N 2 and so on: each successive term is suppressed by an extra factor of 1/N C . If, as 
we have assumed before, A (k) does not vanish at k — 0, how can x(k) vanish at k — 0, 
since the extra terms in x{k) are down by powers of l/iV c ? We can answer this question 
by writing x{k) i n a slightly different (yet more explicit) way. 

As any other two-point function of gauge invariant operators, x(k) can be written to 
the lowest order in 1/N C as a sum over one-hadron poles, i.e., one-hadron intermediate 
states (see Ref. @]): 

glueballs n mesons ™ " b n 

where M n is the mass of the n— glueball state and m n the mass of the n— meson. In this 
formula a n is the matrix element of Q(0) to create the n— glueball state: 



a n = (0\Q(0)\n- glueball), (4.8) 

and c n /\/W c is the matrix element of Q(0) to create the n— meson: 

Cn /JW c = (0|Q(0)|n & meson). (4.9) 

From the above-mentioned l/iV c -expansion counting rules, we know that a n and c n are 
both of order one. The first term in ( (4.7|) (the sum over glueball states) is Ao(k); the 
second term (the sum over mesons) is the contribution of quark loops. 

We have to explain how the second term, which is suppressed by a factor 1/N C , can 
cancel the first one. Although this cancellation is impossible for a general value of k, it 
can happen at k = (as we require!) if there exists a meson state with a squared mass of 
order 1/N C : since this meson couples to Q(x), it must be a pseudoscalar flavour singlet. 
We call it the rj', since the rj' is the lightest flavour-singlet pseudoscalar in nature. At 
this point, it is also easy to derive a formula for the mass of the rj'. The 7/-term is the 
only meson term in ( |4.7| ) comparable in magnitude to A o (0) (both of them being of order 
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1), and must by itself cancel Aq(0), in order that the total %(0) vanishes. Thus we obtain 
that: 



C '' " 4,(0). (4.10) 



We recall from Eq. flO|) that: 



(0|Q(0)|i/). (4.11) 



We can now use the anomaly equation for the U(l) axial current = J2i=i Qilnlsli (L 
being the number of massless flavours; moreover: 75 = — 27°7 1 7 2 7 3 ): 

&*J$ = 2LQ, (4.12) 

and substitute Q(0) in Eq. (|CT| ): 

CV ' (4.13) 



N c 2L 



We have that (0\d^J^\rj') = -ip»{0\ J$\rf), and by definition {0\J$\rf) = is/Wp^ 
(with y/2L taken out in this way, F v > is independent of the number of flavours, to the 
lowest order in 1/N C ). Combining these formulae, we obtain: (Old 11 J^\rj') = ^/2LmiiF rj i . 
F v i is not known experimentally, but to the lowest order in 1/N C it is equal to the pion 
decay constant: F v r ~ F^ ~ 93 MeV. From all of this we derive that: 

2 -1 



Substituting this into Eq. ( f4.10|) , we finally obtain the formula for the squared mass of 
the rf: 



2LA 



m v> = ~Ev-> ( 4 - 15 ) 

7T 

where A = Aq(0) is just the topological susceptibility in a QCD without quarks, i.e., in a 
pure Yang-Mills theory. 

Let us make some comments about formula ( [4.15| ). From the l/iV c -expansion theory, it 
is known that F n is of order 0(y/Nc) (see Ref. [fflfl : in fact the current-current correlation 
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function, which is of order N c , has a term F 2 /k 2 ). Since A = Aq(0) is of order 0(1) in N c , 
from Eq. (|4.15|) we derive that mi, is of order 0(1/N C ). This just means that the mass of 



the rf vanishes in the limit N c — > oo. In the same limit also the anomaly vanishes, since 
we have that: 

The fact that the anomaly vanishes when N c —>■ oo is perfectly natural, since the anomaly 
comes from a diagram with one quark loop, and quark loops are of order 1/N C . At 
N c = oo, the U(l) current is conserved and (unless there is a physical axial symmetry) 
there must be a massless U{\) Goldstone boson: it is the rf. At finite (but large) N c , it 
receives a mass of order 1/N C : also this seems to be natural, because the squared mass 
of an approximate Goldstone boson is linear in the symmetry breaking parameter. The 
quantity A has been determined by Monte Carlo simulations in lattice pure-gauge QCD, 
obtaining a result of the right order of magnitude (see, e.g., Refs. ]48], (49| and [|15|j): 



A ~ (150 + 200 MeV) 4 . (4.17) 



5. Veneziano's mechanism 



Just after Witten's proposal that the U(l) problem could be solved in l/iV c -expanded 
QCD, G. Veneziano showed |10[ that this picture comes out to be automatically consistent 
with expected ^-dependences and anomalous Ward Identities, if a (modified) Kogut- 
Susskind mechanism is used. This is what we call "Veneziano's mechanism", which we 
now describe in its main aspects. 

We must first formulate precisely the U(l) problem in 1/N C language and discuss 
which is the relevant large-iV c limit to be taken into account. Let us consider a QCD 
Lagrangian with an SU(N C ) colour gauge group and Nf quark flavours. L of these are 
taken to be light quarks with common mass m <ti Aqcd, while the remaining ones have 
masses greater than Aq C d (Aqcd is, as usual, the renormalization group invariant scale 
of QCD: in the MS renormalization scheme, for example, Aqcd ~ 0.5 GeV). 
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We consider the Topological Expansion (TE) of QCD discussed in Ref. |29[ and ob- 
tained by taking N c , Nf — > oo with a fixed ratio and with fixed g 2 N c . If at the leading 
(planar) level we have both confinement and the Nambu-Goldstone (NG) mechanism, 
then there is an L 2 -plet of light NG pseudoscalars of mass: 

m NS = F Ns(- 2m Qaqa) = i^|(A 2 ), 

a—l,2,...,L, no sum over a, (5.1) 



2 

NS 



where NS stands for "non-singlef of SU (L) flavour] F^s is the analogous of F T ; both F J 
and (qq) are of order N c |28|, |29], |32|, so that m 2 NS oc ttiKqcd- At this planar level m^s 
is also the mass of the SU(L) singlet, because there are no OZI violating qq annihilation 
diagrams at leading order in 1/N C : the L 2 physical mesons are diagonal in the qq basis 
and all degenerate. Another consequence of planarity is that decay couplings (F n , F K , . . .) 
are all related. We find that: 

{NS— pseudoscalar, a|A (0)|0) = —iS a bq^F NS , 
(5-pseudoscalar | J^(0)|0) = -iq^V2LF NS , (5.2) 

where the = g7 At 7s|A a g are the SU(L) axial currents and jQ = J2i=i qil^l^qi is the 
U(l) axial current. All this happens at the planar level. When going to the next-to- 
leading (non planar) diagrams, we must find that the L 2 degeneracy is broken down to 
an L 2 — 1 degeneracy of NG bosons together with a heavier SU(L) singlet of finite mass, 
even in the chiral limit m (^yvs) — * 0. Since we are taking the limit of light quark masses 
(itl/Aqcd —* 0), it turns out that, to have a smooth limit, it is better to consider the TE 
expansion of QCD, rather than the simple N c — > oo limit of 'tHooft. In this TE expansion 
one keeps L/N c fixed with inclusion of an arbitrary number of quark loops at each order, 
and then one lets (for simplicity) L/N c — > 0. More precisely, the limit that one considers 
is the following: 

' 2 !S ~ mA QCD < J^ A QCD < A QCD- ( 5 -3) 



NS ""^QCV ^ Jj-^QCD ^ IV QCD- 



In this limit one can saturate algebrically the Ward Identities (WI's) of the theory. The 
relevant WI's of QCD have the form [|33| : 

J d 4 x(Td»A;(x)d»A b u (0)) = -i5 ab {A 2 ), (5.4) 
2L J d A x(TQ(x)D (L) (0)) = -2tL(A 2 ) - J d 4 x(TD {L) (x)D {L) (0)), (5.5) 
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AL 2 J d A x(TQ(x)Q(0)) = 2iL(A 2 ) + J d A x(TD {L) (x)D {L) (0)). (5.6) 

Here = 2im Y%=\ QiJsli is the mass term appearing in the divergence of the U(l) axial 
current: d^J^fhx) = D(l)(x) + 2LQ(x); and A 2 = —2mq a q a (no sum over a = 1, 2, ... , L). 
Eq. (|5.4j) is the usual anomaly-free WI for the .SL^L) -non-singlet currents: it implies the 
existence of NG bosons of mass m N s given by ( |5.1| ). 

We will concentrate on the two anomalous WI's, ( |5.5|) and ( |5.6|) . Let us assume, as 
proposed by Witten (see Section 4), that the left-hand side of Eq. ( |5.6| ) is non-zero at 
the leading order in 1/N C of pure Yang-Mills (YM) theory. To obtain this, we introduce 
an axial four-vector ghost of propagator —ig^v/q 2 , coupled to the current of Eq. ( |3.2|) 
with strength iX Y M(q 2 ), where \ Y m = 0(1). Therefore, in pure YM and large N c , the 
left-hand side of Eq. (|5.6j ) may be written as: 



AL 2 {i\ YM {q 2 )f{iq^){-iq v )^ \ q ^ = AtL 2 X 2 YM (0) = 0(L 2 ). (5.7) 



Q 2 



Instead, the right-hand side is of order 0(mLN c ) (in the absence of SU(L) singlets of 
squared mass 0{m)). To eliminate this discrepancy, we have to make a resummation 
of quark loops: in fact, even if they are suppressed by factors of (L/N c ) n , they con- 
tain intermediate states of squared mass of order 0(m). Thus we also assume that the 
ghost couples to the (planar) SU(L) singlet (which we will denote by t]l) with strength 
i^L/N c X v (q 2 )q^ (with again independent of L, N c and m). Then we can sum all quark 
loops in the leading non-planar (cylinder) topology, obtaining the following expression 
for the right-hand side (RHS) of Eq. (E] 



AiL X YM (q 



- m Ns) | 

Q 2 ~ m 2 NS - ±X 2 (q 2 ) m 2 NS + ^(0) 



NS 



Ai L 2 X 2 YM (q 2 )(q 2 -m 2 NS ) AiL 2 X 2 YM (0)m 2 NS 



As a consequence of the resummation, now the RHS of Eq. ( |5.6|) is of order 0(mLN c ), 
as required by the WI (remember that m 2 NS = 0(m)). When considering the two-point 
function of Q(x) at four-momentum q, we find that: 



J q 1 - (mira + -4r Xi (a 2 )) 
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Using the expansion (|4.7| ) given in Section 4, one immediately derives from Eq. Q5.9Q that 
the r/^-mass 777.5 is given by: 



m \ L = m s = m NS + j^ x2 r 1 ( m2 s)- ( 5 - 10 ) 
Furthermore, by going on the r/^-pole and using Eq. (|5.2|) , one obtains that: 

(0\Q(0)\vl) = ^=|m s = X Y M(ml)X v (ml)J^-. (5.11) 



Combining Eqs. (|5.10|) and Q5.11|) , we get, in the chiral limit m {m 2 NS ) — > 0: 



m s = 2L \^ m s> . (5.12) 



We can assume that the g-dependences of Xym and are really negligible, since the range 
of q 2 to be considered is < q 2 < (L/N c )Aq CD <C Aq CD an d variations in this range are 
non-leading in L/N c , if we assume no massless glueballs. (A lattice confirmation of the 
g-independence of Xym was given in Refs. and [0.) Therefore we may substitute 
everywhere Xym{^s) anc ^ \{ m %) with their values calculated at q 2 = 0: Xym(0) and 
A^(0). In this way, Eq. (|5.12| ) acquires exactly the form of Eq. ( |4.15|) derived by Witten: 



m 2 s = 2 J^M, (5.13) 

*NS 



Ay M (0) being precisely the YM topological susceptibility (see Eq. ( |5.7|) ). Using Eqs. ( |5.8|) , 
( |5.10|) and ( |5.13| ), one can easily verify that the WI (|5.6|) is saturated (i.e., its left-hand 



side and its right-hand side turn out to be equal already at this level of approximation): 

ULN c m% s (^^) =2iLm 2 NS F 2 NS . (5.14) 

And the same happens also for the other two WI's. In the above-described picture we 
have considered an ideal world where we have N c colours and L light flavours of common 
mass m <C Aqcd- In the real world we have iV c = 3 colours and L = 3 light flavours, with 
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different masses: m u , ma, m s <C Aqcd (for m.s this relation is less strong that for m u 
and rrid). Nevertheless one can use the above results (in particular Eq. ( 5.10|) ) to derive 



the squared mass matrix for the pseudoscalar nonet (ir, K,rj,rj'). What one finds at the 
end is the following very interesting formula for the ?/-mass: 

6A 

— = m\, + m 2 v - 2m 2 K , (5.15) 

where A = Ay M (0) is the YM topological susceptibility. Inserting here the experimental 
values of the meson masses and of the pion decay constant (F n ~ 93 MeV), one derives 
that A is of order A ~ (180 MeV) 4 . Lattice calculations have confirmed such an order of 
magnitude for A (see, e.g., Refs. ER M and fllSf). 



6. The phase transitions of QCD 



All the above refers to the theory at T = 0. In the following of the paper, we discuss the 
role of the U(l) axial symmetry for the phase structure of QCD at finite temperature. 
One expects that, above a certain critical temperature, also the U(l) axial symmetry will 
be restored. We will try to see if this transition has (or has not) anything to do with the 
usual chiral transition. Let us define the following temperatures: 

• T c h: the temperature at which the chiral condensate (qq) goes to zero. The chiral 
symmetry SU(L) eg) SU (L) is spontaneously broken below T c h and it is restored 
above T c h. 

• Tuny, the temperature at which the U(l) axial symmetry is (approximately) re- 
stored. If (qq) also the U(l) axial symmetry is broken, i.e., the chiral condensate 
is an order parameter also for the U(l) axial symmetry. Therefore we must have: 
Tu(i) > T ch . 

• T x : the temperature at which the pure-gauge topological susceptibility A (approxi- 
mately) drops to zero. Present lattice results indicate that T x > T ch |15|]. Moreover, 
the Witten-Veneziano mechanism implies that Tym > T x . 
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The following scenario, that we will call " SCENARIO 1 ", in which T x < T c h, is, therefore, 
immediately ruled out. In this case, in the range of temperatures between T x and T ch 
the U(l) axial symmetry is still broken by the chiral condensate, but the anomaly effects 
are absent. In other words, in this range of temperatures the U(l) axial symmetry is 
spontaneously broken (d la Goldstone) and the rf is the corresponding NG boson, i.e., it is 
massless in the chiral limit sup(mj) — > 0, or, at least, as light as the pion tt, when including 
the quark masses. This scenario was first discussed (and indeed really supported!) in Ref. 
37| . It is known that the U(l) chiral anomaly effects are related with instantons ||. It is 



also known that at high temperature T the instanton-induced amplitudes are suppressed 



by the so-called "Pisarski-Yaffe suppression factor" |38|, due to the Debye-type screening: 



dMinstiT) ~ dMinstiT = 0) • exp 



9 9 , /2N C + L 
-vrVT 2 [S- 



(6.1) 



p being the instanton radius. The argument of Pisarski and Wilczek in Ref. |37| was t ne 
following: "If instantons themselves are the primary chiral-symmetry-breaking mecha- 
nism, then it is very difficult to imagine the unsuppressed U{1)a amplitude at T ch ." So, 
what was wrong in their argument? The problem is that Eq. ( |6.1|) can be applied only 
in the quark-gluon plasma phase, since the Debye screening is absent below T ch . Indeed, 
Eq. ( |6.1| ) is applicable only for T > 2T ch and one finds instanton suppression by at least 
two orders of magnitude at T ~ 2T c h (see Ref. [17| and references therein). Moreover, the 
qualitative picture of instanton-driven chiral symmetry restoration which is nowadays 
accepted has significantly changed since the days of Ref. [|37|]. It is now believed (see 
Ref. [0 and references therein) that the suppression of instantons is not the only way 
to "kill" the quark condensate. Not only the number of instantons is important, but also 
their relative positions and orientations. Above T c h, instantons and anti- instantons can 
be rearranged into some finite clusters with zero topological charge, such as well-formed 
"instanton-anti-instanton molecules" . 

Therefore, we are left essentially with the two following scenarios. 



SCENARIO 2 : T ch < T u{1) , with T, 



eh 



T ~ 
± x 



T, 



[/(!)• 



If T, 



eh 



T 

- 1 V 



then, in the case 



of L = 2 light flavours, the restored symmetry across the transition is U(1)a ® SU(2)l <g> 
SU(2) R ~ 0(2) ® 0(4), which may yield a first-order phase transition (see, for example, 
Ref. 0). 

SCENARIO 3 : T c h that is, the complete U(L)i ® U(L)r chiral symmetry is 

restored only well inside the quark-gluon plasma domain. In the case of L = 2 light 
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flavours, we then have at T = T c h the restoration of SU (2)l ® SU (2) R ~ 0(4). Therefore, 
we can have a second-order phase transition with the 0(4) critical exponents. L = 2 
QCD at T ~ T ch and the 0(4) spin system should belong to the same universality class. 
An effective Lagrangian describing the softest modes is essentially the Gell-Mann-Levy 
linear sigma model, the same as for the 0(4) spin systems (see Ref. [0). If this scenario 
is true, one should find the 0(4) critical indices for the quark condensate and the specific 
heat: (qq) ~ \{T - T ch )/T ch \ - 38±om , and C(T) ~ |(T - T ch )/T ch \ OA9±om . Present lattice 
data partially support these results. 



7. The U(l) chiral order parameter 



We make the assumption (discussed in the Introduction and in the previous Section) that 
the U(l) chiral symmetry is broken independently from the SU(L) £g> SU(L) symmetry. 
The usual chiral order parameter (qq) is an order parameter both for SU(L) ® SU (L) 
and for U(1)a- when it is different from zero, SU(L) ® SU(L) is broken down to SU(L)y 
("V" stands for "vectorial") and also U(1)a is broken. Thus we need another quantity 



which could be an order parameter only for the U(l) chiral symmetry [23|, 40|. 

The most simple quantity of this kind was found by 'tHooft in 1976 (see ||) studying the 
effective Lagrangian of instantons (so that, for historical reasons, we will call this quantity 
"£ e //"). For a theory with L light quark flavours (of mass rrii <^ Aq C d', i = 1, ■ ■ ■ , L), it 
is a 2L-fermion interaction that has the chiral transformation properties of: 



4g - det 



1 + 75 



+ h.c. = det (q sR qt L ) + det {q sL qtR) , (7.1) 

st st 



where s,t = 1, . . .,L are flavour indices, but the colour indices are arranged in a more 
general way (see below); as usual, qi = |(1 + 7s)g and qn = |(1 — 7s)g, with 75 = 
— i7°7 1 7 2 7 3 , denote respectively the left-handed and the right-handed quark fields. Since 
under chiral U(L) ® U (L) transformations the quark fields transform as follows: 

<?L -> V L q L , q R -> V R q R , (7.2) 



21 



where Vl and Vr are arbitrary L x L unitary matrices, we immediately derive the trans- 
formation property of 4// under U(L) <g> U(L): 

U(L) <g> U(L) : 4// -> det (^) det(VR)*det (q sR q tL ) + h.c. (7.3) 

This just means that 4// is invariant under SU(L) cg> SU(L) ® Z7(l)y, while it is not 
invariant under U(1)a'- 

U(1) A : £g} - e~ i2La det (q sR q tL ) + h.c., (7.4) 

where a is the parameter of the U(l) chiral transformation (U(I)a'- q — » e~* a75 g, i.e., 
<?l -> e~ J °g L and g fl -> e ,a g fl ). 

As an example let us consider the most simple case, that is L = 2, but with a general 
colour group SU(N C ). It is not hard to find (using the Fierz relations both for the spinorial 
matrices and the SU(N C ) generators in their fundamental representation) that the most 
general colour-singlet, hermitian and P-invariant local quantity (without derivatives) 
which has the required chiral transformation properties is just the following four-fermion 
local operator: 

4^7 2) (a , A)) = F%(a , (3 )e st (<f lR q h sL ■ ? m qt L + <? 1L q h sR ■ q c 2L qf R ) , (7.5) 
where the "colour tensor" F^ c (a , A)) * s gi ven by: 

F%(a , (3 ) = a « + A)«, (7.6) 



a and /3 being arbitrary real parameters. In Eq. (|7.5|), a,b,c,d = 1, . . . , N c are colour 
indices; s and t are flavour indices and e st = —e ts , e 12 = 1. Dirac indices are contracted 
between the first and the second fermion field and also between the third and the fourth 
one. Note that if we choose ao = N c and (3q = — 1, C e ^ (ao. A)) just becomes (up to a 
proportionality constant) the effective Lagrangian found by 'tHooft in |5| to describe the 
instantons. 



the vacuum expectation value of C ef J («o, A> 



Now, to obtain an order parameter for the U(l) chiral symmetry, one can simply take 

(L=! 
eff 

4« 2) («o, A)) = (4// 2) ("o,/9o))- (7.7) 
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The arbitrarity in the choice of ao and (3q (indeed of only one of them, since only their 
ratio is relevant) can be removed if we require that the new U(l) chiral condensate is 
"independent", in a sense which will be explained below, of the usual chiral condensate 
(qq)- As it was pointed out by Shifman, Vainshtein and Zakharov in PTJ, a matrix element 



of the form (qTiq ■ qT 2 q) has a big contribution proportional to the square of the vacuum 
expectation value (v.e.v.) of qq. This contribution corresponds to retaining the vacuum 
intermediate state in all the channels and neglecting the contributions of all the other 
states; we call this contribution the "disconnected part" of the original matrix element: 

(qT iq ■ qT 2 q) dtsc . = ± [(TrT 1 ■ TrT 2 ) - TrQ^)] (qq) 2 , (7.8) 

where the normalization factor N is defined as: 

(QaQb) = (7.9) 

(that is N = Saa] qq = J2a QaQa) an d the subscripts A, B include spin, colour and flavour. 
When considering the operator C e fJ (ao, /?o) defined in Eqs. (|7.5|) and ( |7.6| ), we find the 
following expression for its disconnected part: 

(4/? 2) K A>)W = ^[N c (2a + fa) + (a + 2{3 )](qq)\ (7.10) 

(where: (qq) = (uu) + (dd)). From this last equation we immediately see that the discon- 
nected part of the condensate CLjj (cyq, Po) vanishes with the following particular choice 
of the coefficients cto and (3q (only their ratio is really relevant): 

^ = -^4- (7.11) 

a N c + 2 1 ; 

In other words, the condensate Eq. ( |7.7| ) with ao and (3q satisfying the constraint Eq. 
( 7. 11|) does not take contributions from the usual chiral condensate (qq). To summarize, a 



good choice for a (7(1) chiral condensate which is really "independent" of the usual chiral 
condensate (qq) is the following one (apart from an irrelevant multiplicative constant): 

C^= 2) = ((m - ^±l8 a d 8^e st (q a m q b sL ■ q c 2R qf L + q a 1L q b sR ■ q c 2L q d tR )). (7.12) 
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As a remark, we observe that the condensate C^-m so defined turns out to be of order 
0(g 2 N 2 ) = 0(N C ) in the large- N c expansion (this will be derived also in Section 8 by 
simply requiring that the 1/N C expansion of the relevant QCD Ward Identities remains 
well defined when including this new condensate |2"2|). In the particular (yet physical!) 
case N c = 3, the condensate Eq. (|7.12|) becomes: 

4w 2) = am - \my l • + ■ (7.13) 

So far we have considered the most simple case L = 2. However, this procedure can be 
easily generalized to every L, and we can take as an order parameter for the U(l) chiral 
symmetry: 

= (£$). (7-14) 

As we have done in the case L = 2, the colour indices may be arranged in such a way that 
the U(l) chiral condensate does not take contributions from the usual chiral condensate 
(qq): as a consenquence of this, the new condensate will be of order 0(g 2L ~ 2 N^) = 0(N C ) 
in the large- N c expansion. The case of the real world is that in which L = 3 (three light 
flavours, u, d and s, with masses m u , rrid, m s small compared to the QCD mass-scale 
Aqcd)- In the next Sections we will try to derive which are the physical consequences 
of assuming that ^ up to a temperature Tv{\) far above the usual chiral critical 

temperature T c ^. For doing this, we will analyse the relevant QCD Ward Identities by 
making an expansion in 1/N C and a chiral expansion in the light quark masses. In Section 
9 we will see how to derive a new chiral effective Lagrangian which also incorporates the 
U(l) chiral order parameter. 



8. Saturation of Ward Identities 



In this Section we make an analysis of the relevant QCD Ward Identities (WI's). Let us 
consider a set of ideal worlds described by the usual QCD Lagrangian with an SU(N C ) 
colour gauge group and Nf quark flavours, where N c and Nf are considered as variables. 
Among these Nf quarks, L are taken to be light quarks with mass mi <C Aq C d, A-qcd 
being the renormalization-group invariant scale of QCD. The saturation of the relevant 
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QCD Ward Identities will be obtained in a 1/N C expansion, N c being the number of 
colours p8[ ^ . 



When deriving the relevant WI's as in Refs. j|2| and [10|, one only uses the canonical 
(anti-) commutation relations at equal times: 

which do not absolutely depend on the temperature T of the physical system: they are a 
fundamental law of quantum mechanics. For this reason the relevant WI's at T 7^ are 
simply derived from those at T = by substituting everywhere the vacuum expectation 
values with the usual thermal average over the Gibbs ensemble (/3 = 1/7 1 ): 

(0|o|0> ^l^W = " ' (Q> ' (8 ' 2) 

n 

H being the total hamiltonian of the system, E n and |n) its eigenvalues and eigenstates 
(H\n) = E n \n)). 

In the case of the SU(L) £g> SU(L) chiral symmetry, one immediately verifies that: 

(iQ a A (x ),tq l5 T b q(y)} Xo=yo ) = -t5 ab ^(qq), (8.3) 

where Q\(xq) is the charge operator for the SU(L) chiral symmetry (i.e., Qa( x o) 
j d 3 Xv4g(x, x ), where A® = q^/y^^q is the SU(L) axial current). From Eq. ( |3.3| ) one 
derives the following WI: 

d i x(Td^A^(x)tq l5 T b q(0)) = t5 ab j(qq), (8.4) 



If (qq) 7^ (in the chiral limit sup(mj) — > 0), the anomaly-free WI ( |8.4| ) implies the 
existence of L 2 — 1 non-singlet NG bosons, interpolated by the hermitian fields O b = 
%q^T h q. Similarly, in the case of the U(l) axial symmetry, one finds that: 

([Q5 ( x o),iqj5q(y)]x =y ) = _2i (^)> ( 8 - 5 ) 
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where Q 5 (xo) is the charge operator for the U(l) chiral symmetry (i.e., Q 5 (xo) = 

x ), where = q^f/y 5 q is the usual U(l) axial current). From Eq. 
one derives the following WI: 

J d 4 x(Td»J^(x)iq-f 5 q(0)) = 2i(qq). (8.6) 

But this is not the whole story! One also derives that: 

[Q^\x ), O^\y)} X0=yQ = -2H • C$ f {y), (8.7) 

where ~ det(q S Rq t L) + det(q S Lq t n) is the 2L-fermion operator discussed in the pre- 
vious Section and the hermitian field is so defined: 

OP ~ i[det(q sR q tL ) - det{q sL q tR )]. (8.8) 

From Eq. (|8.7D one derives the following WI: 

/ d'x(Td^(x)O P L \0)) = 2Li(4g(0)>. (8.9) 

Remembering that d^J^) = + 2LQ, where D^l) = ^Y,i=i m iqn^qi an d Q(x) is the 
topological charge density, Eq. ( |3.9|) can also be written in this way: 

2L J d 4 x(TQ{x)O P L) {0)) = 2Li(4//(0)) - / ^(TD^^O^^O)}. (8.10) 

If the U(l) chiral symmetry is still broken above T c h, i.e., (£^(0)) ^ for T > T c h (while 
{qq) — for T > T c h), then this WI implies the existence of a (pseudo-)Goldstone boson 
(in the large- N c limit!) coming from this breaking and interpolated by the hermitian field 
Op. Therefore, the U{1)a (pseudo-)NG boson (i.e., the ?/) is an "exotic" 2L-fermion state 
for T > T c /j! The fact that (£^(0)) is different from zero above T c ^ in the chiral limit 
of zero quark masses, forces us to introduce a new 2L-fermion effective vertex with a 
coupling constant that we may call "7". This is shown schematically in Fig. 1. 

Now let us consider the implications of this picture on the saturation procedure of 
the other relevant QCD Ward Identities in the region of temperatures between T c h and 
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Figure 1: The new 2L-fermion effective vertex. 

T x < having assumed that T ch <T X < W^ - 111 other words, we suppose to be 

in a region of temperatures T > T c h where the SU(L)®SU(L) chiral symmetry is restored 
(i.e., (qiqi) = 0, in the chiral limit sup(m/) — > 0), the topological susceptibility of the pure 
Yang-Mills (YM) theory is different from zero (i.e., (QQ)\ym = ^ym 0i nere an d i n 
the following we make use of the notation: (AB) = j d 4 x(TA(x)B(0))) and the U(l) 
axial symmetry is still broken. Following Veneziano's model for the resolution of the U(l) 
problem [JITJ, we may get that (QQ)\ym = iAy M is different from zero by introducing an 
axial four-vector ghost of propagator —ig^v/q 2 coupled with strength i\ym to the current 
Kp, whose divergence is the topological charge density Q (Q = d^K^): 



(QQ)\ym = {(i^ym) {iq^) — 5— } =^ym- ( 8 - n 




In order to explain, in a 1/N C picture (see Ref. 0), how it can be that (QQ) = 0, in the 
chiral limit sup(m z ) — > 0, while (QQ)\ym is of order 0(1) 7^ with respect to the light 
quark masses, we have to consider (see Ref. ]ITJ) that there are intermediate pseudo- 
scalar bound states which couple to the ghost and which therefore must be resummed. 
Since we are above the SU(L) ® SU(L) chiral phase transition, it is natural to suppose 
that there is just one S'L r (L)-singlet pseudo-scalar state which plays such a role and which 
also breaks the U(l) chiral symmetry. Following Section 7, we take as an order parameter 
for the U(l) chiral symmetry alone the expectation value of the 2L-fermion interaction 
term of the type of the one introduced by 'tHooft in Ref. 0. 

With the notation introduced before, the relevant WI's at T 7^ take the form: 



d A x{TDi(x)Q{Q)) = -2 / d 4 x(TQ(x)Q(0)), (8.12) 
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4 J d 4 x(TQ(x)Q(0)) = 2%(A,(0)) + / d 4 x(TA(x)L>,(0)). (8.13) 

where: A;(x) = —2miq~iqi and A(^) = 2%miqi^q\ (so that D(n = Z^i-Dj). We require 
their saturation at the leading order in a l/iV c expansion. To have a smooth limit when 
iV c — > +00, it is necessary not only to take the QCD coupling constant "g" of the form 
9 = 9o/ V^Vc with g of order 0(1) with respect to iV c (see Refs. [^8|] and [p9| ) but also to 
take the new coupling constant "7" of the form: 

7 = j^k, (8-14) 

with 70 still of order 0(1) in N c . One can easily be convinced of this by analysing the in- 
dependence of the relevant Feynman diagrams for Eqs. (|8.12 ) and (|8.13|) , considering also 



all possible insertions of the new effective vertex "7" . On the other hand, Eq. ( |S.14 ), with 
7o of order 0(1) in N c , is guaranteed by a proper choice of the operator Cjff, such that its 
expectation value C^l-, = (£e//) comes out to be really independent of the usual chiral 
condensate (qq), as discussed at the end of Section 7: in fact, this implies C^)^ = 0(N C ), 
from which one immediately derives Eq. ( |8.14j ), after a comparison with Fig. 1. 



2im,7 




Q + 2im,7 




Figure 2: The LHS of Eq. (fOg ). 



Let us consider the WI ( ^.12| ), which we briefly write as (DiQ) = —2(QQ). At the 
leading order in 1/N C the left-hand side (LHS) term is schematically of the form given in 
Fig. 2. Shaded blobs represent the sum of all diagrams of a given topology (i.e., what we 
may call a ll gluonic covering" of the basic diagrams). The first piece is of order 0(1) in 
N c and of order 0(mf) with respect to the light quark masses; the second one is of order 
0{p(N^~ x ) = 0(1) in N c and of order 0(11^=1 m k) with respect to the quark masses. The 
right-hand side (RHS) of Eq. ( |3.12j ) at the leading order in 1/N C is represented in Fig. 3. 
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Figure 3: The RHS of Eq. fl812l) . 

It is of order 0(1) in N c (i.e., X Y m = 0(1)) and of order 0(1) with respect to the quark 
masses. How can one match the (LHS) with the (RHS) of the WI (|8.12|) at the leading 



order in 1/N C 7 We observe that the second diagram in Fig. 2 contains, as an intermediate 
state, just the pseudoscalar SU (L)-smglet boson described by the field <y p\x). This is a 
new "exotic" particle which has essentially the same quantum numbers (spin, parity and 
so on) as the rj at T = 0: yet it has a completely different quark content, being a 2L- 
fermion field of the form ~ i[det(q sL q tR ) — &et(q S R<itL)]- We may obtain that the WI ( |8.12| ) 
is saturated by this bound state \Sx), assuming that it couples to the four-vector ghost 
of Kn with strength if s }- — q^ (g M being the four-momentum), where f s is of order Oil) 

in N c . (The fact that f s is of order 0(1) in N c will be verified below, as a consequence of 
Eq. ( |8.16p .) It also couples to the operator Di with strength m^F Dl (observe that Di is 
the divergence of a four-current: D\ = d^(J l b4l — 2K^), where J l 541 = qiJ^Jsqi)- 

(Q\Di\S x ) =m 2 F Dl , (8.15) 

where mo is the mass of the state \Sx) at this OZI order of the expansion (i.e., before 
doing the resummation). We obtain that: 

/ / 1 \ ///'"' 



LHS(8.12) = { (F Dl m 2 )- ^ if*^=% -^^ymHQu) 



-iFoi — 7== Ayj^ . (8.16) 



Fd 1 is of order O(y'yN^) = 0(yN c ). This can be seen considering the relevant Feynman 
diagrams, in the correlation function (DiDi), which contain the state \Sx) as an inter- 
mediate (propagating) state (the procedure is the same which is also used in the usual 
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" ' quark-anti-quarW meson case, and it is well described in the first part of Ref. 
these relevant diagrams contain two 7-vertex insertions and, by virtue of Eq. fl8.14j) , they 
are of order 0(N C ) with respect to N c . Since, as we have seen above, Eq. ( |8.16| ) is of 
order 0(1) with respect to N c , we immediately derive that f s is of order 0(1). 

Comparing the expression ( |3.16D with the RHS of Eq. ( |8.12| ), written in Fig. 3, we 
derive that: 



D, 



fs 



(8.17) 



As it must be, if we want to saturate the WI (|8.12|) , the coupling constant has the 
same value (Fjj) for every index I — 1, 2, . . . , L. From the structure of the second piece of 
Fig. 2 and from Eq. ( |8.15| ), one immediately recognizes that the squared OZI mass m\ is 
proportional to the product of the L light quark masses (assuming, as usual, that F Dl is 
of order 0(1) with respect to them): 



m 



All of this comes directly from the leading-order terms of the expansion. All quark loops, 
inserted in the leading non-OZI topology, can then be easily resummed as a geometric 
series and we get that: 

{QQ) KT iq) = I dV x (TQ(x)Q(0)) = 



i\ 2 YM (1 + {* - 1 2 + ...)= i\ 2 YM q - m ° f2 . (8.19) 



Also, in the same way: 

{DiQ) F ' T '(q) = I d 4 xe^(TA(x)Q(0)) 



Dll7l ° y^iVj YM( 1 2 -»o\ q 2 -ml " ')" 

fs 1 

= iF Dl ml-j==\ YM — — — h— . (8.20) 

The right-hand sides of Eqs. ( |8.19|) and ( |8.20| ) have a pole corresponding to the physical 
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mass of the pseudoscalar singlet boson |<Sx): 



f 2 A\ 2 
rnl = m l + ^ = rnl + ^. (8.21) 



(In the last passage we have used Eq. (|8.17| ).) So the boson \Sx) acquires a squared 



"topological" mass fs/jN^ = ^ymZ-^dj which is different from zero also in the chiral 
limit of zero quark masses (when, on the contrary, the OZI mass mo vanishes). Moreover 
this squared mass is of order 0(1/N C ) in the number of colours, just like the anomaly 



term and the rj squared mass in the Witten-Veneziano model |9|, |10|: in this sense we 
may say that this new "exotic" state is a ll lighf state, at large N c . This result will also 
be obtained in Section 14 in an effective Lagrangian model. It is also easy to verify, using 



Eq. (ETT7D , that the expressions (|87T9D and flOO)) , in the limit q -> 0, satisfy the WI 



(ETI3) (that is: (D t Q) = -2(QQ)); they are both of order O^N^) = 0(N C ) with respect 



to N c and of order 0(Y[^ =1 m k ) with respect to the light quark masses. In fact, using Eqs. 
(1817D and (gl8|) , we find: 



(QQ) = (QQ) RT iq - 0) = i\ 2 YM 



— > '^M^Tfmo = ^— ™n = 01 I I m fc I . (8.22) 

sup(m ; )^0 J s 

Also this result is in total agreement with what we will find in Section 14 using an effective 
Lagrangian model. And what about the chiral condensate (<&<&)? At the leading order 
in our expansion in N c and in the light quark masses mi (I = 1, 2, . . . , L) it receives two 
different contributions of order 0(N C ): one from a mass insertion and another from a 
7-vertex insertion (see Fig. 4). 

Also this last behavior of the chiral condensate with respect to the light-quark masses 
will be found in Section 14. (In the Appendix we also derive the same quark-masses 
dependence for the topological susceptibility and the chiral condensate using a simple 
instantonic model.) As a curiosity, the fact that the chiral condensate (qiqi) does not 
vanish as the mass mi goes to zero (since there remains the piece O2 proportional to the 
product of the remaining L—l quark masses) could appear to be rather unnatural. Yet it is 
not so strange, if one considers that the expectation value (Cj^ JO)) is an order parameter 
also of the C/(l)^ chiral symmetry for the flavour (i.e., the group of transformations 
qi — > e~ ian5 qi } which rotates only the flavour q t ). Therefore, if (£^1(0)) 7^ 0, as in our 



L 
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= 0,(171.) + o 2 (n k/l m k ) 



Figure 4: The chiral condensate above T< 



eh ■ 



case, the U(l)^ symmetry is broken even when the mass mi of the flavour q\ goes to zero. 



Finally, we can immediately verify that also the second relevant WI Q8.13Q is saturated 



4<QQ> = 2i<%(A(0)) + (DiDj). (8.23) 

Let us take I = j in the last equation. Then it is easy to recognize that the 0\(mi) term 
in Fig. 4, that is the part of order 0(N C ) in iV c and of order 0(mf) in the quark masses 
for (Aj(0)), cancels exactly with the analogous part of (DiDi): 

2t(M0))\N c , m f + (AA>k, m? = 0. (8.24) 

The part of (Aj(0)) of order 0(^N^) = 0(N C ) in N c and of order 0(J\k=i m k ) in the quark 
masses, coming from the C^dlfc^ m k) term in Fig. 4, is forced to be directly proportional 
to the "full" topological susceptibility (since (DiD\)\ » rL rr L = 0) : 

7" c 'llfc=i mfc 



2z(A,(0)) \ iN ,^ m = -Aim&AXlmA =A(QQ). (8.25) 

That is: (QQ) = — ifTtyC^dlfc^z m k)- Also this relation will be found in Section 14, as a 
consequence of an effective Lagrangian model. For / ^ j the result is the same, since one 
can easily verify that: 



V/ 



*j : 2i(A { (0)> |^, ntimfc = (A^)l 7jv , )nLimfc , (8.26) 
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in agreement with the fact that both of them must be equal to A(QQ). 



9. The new chiral effective Lagrangian 

We will see in this Section how the proposed scenario, in which the U(l) axial symmetry is 
still broken above the chiral transition, can be consistently reproduced using an effect ive- 



Lagrangian model [21, 22, 23 



It is well known that the low-energy dynamics of the pseudoscalar mesons, including 
the effects due to the anomaly and expanding to the first order in the light quark masses, 
can be described, in the large- N c limit, by this kind of effective Lagrangian [[43], [ 
|§ g7|: 
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C(U, U\ Q) = C (U, U ] ) + -^=Tr[MU + M^U j ] 

2y 2 

+±iQ(x)Tr[\og(U) - log(C/t)] + ^Q 2 (x), 



(9.1) 



where £ describes a kind of linear cx-model: 



C (U, C/t) = hbidJJPtf) - \\lTr[(U^U - p w ■ I) 2 ] 



(9.2) 



I is the identity matrix. A few words about the notation used in Eq. ( |9.1| ). M represents 
the quark mass matrix, which enters in the QCD Lagrangian as SCqcd = —QRMqL+h.c, 
that is: 

(mi \ 

m 2 



M 



V 



(9.3) 



m L 



2 

Q( x ) = &t^ £lJlUpa F{lu(, x )Fpa(, x ) 1S ^e topological charge density and A is a real positive 
coefficient, which is essentially the topological susceptibility in the pure Yang-Mills (YM) 
theory (i.e., in the absence of quarks): 



d 4 x{TQ(x)Q(0))\ YM = iA. 



(9.4) 
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The parameter p w , which is a function of the temperature T, controls the behaviour of 
the theory across the chiral transition, being p n (T < T ch ) = ~F% > and p w (T > T ch ) = 
— \B"l < (p n (T = T c h) = 0). The mesonic field is a L x L complex matrix, which, 
in terms of the quark fields q^ (i = 1, 2, . . . , L) is essentially: 

U U ~ Qj ( 1 + n l5 ) 1i = VjR<liL, (9-5) 



up to a multiplicative constant; as usual, q^ = |(1 + 75)? and q^ = |(1 — 75)9, with 
75 = — i7°7 1 7 2 7 3 , denote the left-handed and the right-handed quark fields. Since under 
chiral U(L) <8> transformations the quark fields transform as follows: 

<?l -> Vbgi, , to -> Vjigfl, (9.6) 

where Vj, and Vr are arbitrary L x L unitary matrices, we immediately derive the trans- 
formation property of U under U (L) <g> U (L) : 

u -> (9.7) 

Let us also observe that C (U, W) is invariant under C/(L) <g> U(L), while the Lagrangian 
£(£/, C/t, Q), in the chiral limit sup(mj) — > (z = 1, 2, . . . , L), is invariant under SU(L) <8> 
SU(L) ®C/(l)y ("V" stands for "vectorial") and it has the same transformation property 
of the massless QCD Lagrangian Cq CD under a £7(1) chiral transformation: 

17(1) A : £gU*) - + 2LaQ(x), (9.8) 

where a is the parameter of the U(l) chiral transformation (U(I)a'- q — > e _ia75 g, i.e., qi — > 
e~ ta qh and g# — > e ia qn). The Lagrangian ( |9.ip considers, as a chiral order parameter, only 
the usual chiral condensate (g?): where it is different from zero, both SU(L) <S>SU(L) and 
U(1)a are broken. If we make the assumption that the U(l) chiral symmetry is restored 
at a temperature greater than T ch , we need another order parameter for the U(l) 
chiral symmetry, the form of which has been discussed in the previous Sections. We must 
now define a field variable X, associated with this new condensate, to be inserted in the 
chiral Lagrangian. The translation from the fundamental quark fields to the effect ive- 
Lagrangian meson fields is done as follows. The operators iq^q and qq entering in the 
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WI ( |S.6|) are essentially equal to (up to a multiplicative constant) i(TrU — TrW) and 
TrU + TrW respectively. Similarly, the operators C e ff ~ det^^g^) + det(q S £,q t R) and 
Op ~ i[det(q S Rq t L) — det(q S Lq t R)] entering in the WI fl8.9|) can be put equal to (up to a 
multiplicative constant) X + X' and i(X — X*) respectively, where 

= det (q sR q tL ) (9.9) 

is the new field variable (up to a multiplicative constant), related to the new U(l) chiral 
condensate, which must be inserted in the chiral effective Lagrangian. Clearly X has the 
following transformation properties under U(L) <g> U(L): 

U{L)®U{L): X -> det(V L )det(V R )*X. (9.10) 

That is, X is invariant under SU{L)®SU{L)®U{1) V , while, under U(1) A , X -> e~ i2La X. 

The new Lagrangian, written in terms of the new fields C/y, X and Q, must have the 
following properties: 

• it must be invariant under U(L) (g> U(L), if one neglects the quark masses and the 
effects due to the anomaly; 

• it must transform in the same way as Cqcd m Q9-81) under the group U(1)a (always 
neglecting the quark masses). 

As usual we will consider the effects of the quark masses by making an expansion in 
powers of M and taking only the linear term. Since the effective Lagrangian has to be 
invariant under a transformation of U as in ( fj.7p and of X as in ( |9 . 1 U| ) together with a 
transformation of M of the form (we recall that M represents the quark mass matrix, 
which enters in the QCD Lagrangian as 5Cqq^ = —q R Mq L + h.c): 

u (9.ii) 



X 



det 

St 



1+75 



Qt 



we see that the most simple term linear in M is always of the following form: 

S£mass = ^| Tr t Mf/ + ( 9 - 12 ) 
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At this point it is rather simple to see that the most simple effective Lagrangian, con- 
structed with the fields U, X and Q and satisfying all the properties listed above, is 
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C(U,U\X,X\Q) = ^TridpUdHJ*) + ^Xd"X^+ 
-V(U,U\X,X^) + -iQ^LUxTrilogiU) - log(t^)] + 



+ \iQ{x){l - ^)[log(X) - log(Xt)] + ^Q\x) 
where the potential term V(U, U*,X,X*) has the form: 



(9.13) 



V(U, U\ X, Xt) = ±%Tr[(rfU - p. • I) 2 ] + \\\ (XtX - p x ) 

-^=Tr[MU + M ] U ] ] - -^[detffltf + det(f/ t )X]. 
2y2 2y2 



(9.14) 



All the parameters appearing in the Lagrangian must be considered as functions of the 
physical temperature T. In particular, the parameters p^ and px are responsible for the 
behaviour of the theory respectively across the SU(L) <g>SU(L) and the U(l) chiral phase 
transitions, according to the following table: 





T<T ch 


T c h <T < T(7(!) 


T > T U{1) 


Ptt 


\Fl > 


-\Bl < 


-\Bl < 


Px 


±F1>0 


|F1>0 


-\B\ < 



[Tab.l] 

(That is: p 7r (T c / l ) = and px(Tu(i)) = 0.) The U(l) chiral symmetry remains broken also 
in the region of temperatures T ch < T < Tf/m, where on the contrary the SU(L) <g>SU(L) 
chiral symmetry is restored. The U(l) chiral symmetry is restored above Tym. Following 
the Introduction, we also assume that the topological susceptibility of the pure YM theory, 
A(T), drops to zero at a temperature T x greater than T c ^. We will see later that, for the 
consistency of the model, it must be T x < In the following it will be sometimes 

useful to integrate out the field variable Q(x) in the Lagrangian fl9.13| ), obtaining: 
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-V(U,tf,X,Xi) + _A{o; 1 TV[log(C/) - Zop(C/t)] + 

+ (l-^)[log(X)-log(Xt)]} 2 . (9.15) 



In order to avoid a singular behaviour of the anomaly term in ( |9.15| ) just above T c ^, where 
the "vacuum" expectation value of U tends to zero in the chiral limit sup(mj) — > and 
where A(T) > 0, we will assume that: 

^(T>T ch ) = 0. (9.16) 



10. Mass spectrum of the theory for T c h < T < Tjj^ 



Let us now consider the Lagrangian ( 9.15 ), where the field variable Q(x) has been in- 
tegrated out; we study the mass spectrum of the theory in the region of temperatures 
T c h < T < Tu(\). As it has been stressed in Ref. jlcj], the linear cr-type model con- 
tains redundant scalar fields, which can be eliminated by taking the limit A 2 — > +oo and 
\ 2 X — > +oo in Eq. ( |9.14 ): so we will expand our results not only in powers of the light 



quark masses rrii (i — 1,2, . . . ,L), but also in powers of 1/A 2 and 1/A^. In the region 
of temperatures that we are considering, we have, according to Table 1, p n = — ~B% < 
and px = \F X > 0. So, because of the form of the potential term ( |9.14| ), in this region 



the SU(L) <8) SU(L) chiral symmetry is restored, being (U'U) = 0, in the chiral limit 
sup(mj) — > 0, while the U(l) chiral symmetry is still broken, being (X^X) = \F\. More- 



over, from Eq. (|9.16|) we have that uj\{T > T c h) = 0. If we parametrize the complex fields 



U and X by separating their real and imaginary parts in the following way: 

Uij = a i:j + ibij ; X = e x + i6 2 , (10.1) 

we immediately find that the minimum of the potential is obtained in the "point" P 
where: _ 

1 ' ' " " 10.2) 
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We first neglect the anomaly term in ( 9.15Q , that is the term proportional to the topological 



susceptibility A of the pure YM theory, and then we study the mass spectrum of the theory, 
by diagonalizing the matrix of the second derivatives in the "point" P of the potential 
V(U, WjXjX^) with respect to the various fields ciij, bij, 6\ and 6 2 - We then find that, 
in this limit, there is a field Sx (which is essentially Q 2 — Q 2 \p = Q 2 ) which has a squared 
mass ml proportional to the product of the light quark masses: 

All the other fields have squared masses equal to X X F X or ^X^F^, so that, in the limit 
A 2 — > +00 and X x — > +00, they are redundant fields. If we now consider the full potential 
V(U,U\X,X^) which appears in Eq. ( p!5|) as the sum of V(U, U\ X, X f ) and of the 
anomaly term: 

V(U, U\ X, Xt) = V(U, U\ X, Xt) - ^[log(X) - log(Xt)] 2 , (10.4) 

we obtain that the field S x acquires a "topological" squared mass equal to In con- 
elusion, the squared mass of the field S x is given by: 

7775 = 777 + —5- > —r. (10.5) 

The physical interpretation of this new particle described by the field S x is rather obvious. 
It is nothing but the would-be Goldstone particle coming from the breaking of the U(l) 
chiral symmetry: in fact, neglecting the anomaly, it has zero mass in the chiral limit of zero 
quark masses. Yet it acquires a "topological" squared mass proportional to the topological 
susceptibility A(T) of the pure YM theory, just as the 7/ in the Witten-Veneziano model 
0, [l0| . And actually it has the same quantum numbers (spin, parity, etc.) as the 7/, 
even if it is a sort of "exotic" matter field of the form ~ ^[det^^g^) — det (q S RqtL)]- Its 
existence could be proved perhaps in the near future by heavy-ions experiments. 

Finally we also observe that, in order to avoid a singular behaviour of 777| x as the tem- 
perature T approaches the U(l) chiral transition temperature Tu(i) (where F x vanishes), 
we must require that the "topological" temperature T x , at which A(T) drops to zero, is 
not greater than Tfjm; that is: 

T c h <T X < Tun) ■ (10.6) 
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11. Mass spectrum of the theory for T < T c h 



In the previous Section we have studied the mass spectrum of the theory for T c h < T < 
It is obviously interesting to see also what happens for T < T c h, where both the 
SU(L) ® SU(L) and the U(l) chiral symmetry are broken. In fact, owing to the form of 
the potential ( |9.14j) , we have, at the leading order in uii, 1/A£ and 1/X X : 



(rfU) = ±I%.I ; (X^X) = \f x . (11.1) 

As we have already explained in the previous Section, we can eliminate the redundant 
scalar fields of the linear cr-type model by taking the limit A^ — > +oo and X 2 X — > +oo in 
Eq. ( |9.14| ). In this limit the potential term gives the following constraints: 

^U=\fI-\ ; X^X = \f\. (11.2) 
So we can take the field U of the form: 




with: 



U = \ ~F v exp , (11.3) 



where the matrices are the L(L— 1) non-diagonal generators of U(L) [or, alternatively: 
$ = Xiali 1 ^a^a + where r a (a = 1, . . . , L 2 — 1) are the generators of the algebra of 
SU(L) in the fundamental representation, with normalization Tr(T a ) = 0, Tr(r a r;,) = 5 a b, 
and SV is the S'C/(L)-singlet field]. And also: 




X = \-Fxexp[^-S x ). (H.5) 



Substituting Eqs. ( |11.3|) and ( |11.5| ) into Eq. (|9.15|) (where the field variable Q(x) has 



been already integrated) and taking only the quadratic part of the Lagrangian, we obtain 



39 



the following expression: 



oL^i f i "oL r 



L 
Cl 

2v/2 



,2 i ,,2 



-7T 7T + 



(3) 



V2 



^^-—Sx) + 



\ - 1 T j=l 



*■ i=i 



where: 



// 



2 -^m 



m, : . 



If we also put, for simplicity: 



ci_ (FA (F^ 
y/2 \V2J \V2, 



2c 



/ 2A(l-wi) 2 +2c 2Ai>i(l-o;i)- 
2Aji(1-cji)-2c 2 _i_ 2Aq^+2c 



2Aji(1-o>i)-2c 



2Awf+2c 
~~Fl 



2Aji(1-oii)-2c \ 
2Ao^+2c 

F2 



2 , 2A^+2c 



(11.6) 
(11.7) 

(11.8) 



we have the following squared mass matrix for the system of fields (Sx, 7Ti, ^2 ,• • • , ttl) : 



11.9) 



It is not hard to find out the L + l eigenvalues of this matrix. In the chiral limit, in which 
the masses m t of the L light quarks are put equal to zero (so that // 2 = 0; i — 1, 2, . . . , L) 
the matrix A has L — 1 null eigenvalues (A = Ai = . . . = A^_ 2 = 0) and two other 
eigenvalues given by: 

A L _i \ Z L t ygj ~ 4Q L 
Al J 2 
where Z L and are so defined: 



(11.10) 



2A[F 2 (1 - coO 2 + LF x ujj] + 2c(F 2 + LF X ) 



Ql = 



ALAc 



(11.11) 
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In other words, we have L — 1 zero-mass states from the ensemble (Sx, ffi, 7T2, • • • , kl) 
which, together with the L(L — 1) zero-mass states n al3 , constitute the L 2 — 1 Goldstone 
bosons coming from the breaking of the SU(L)®SU (L) chiral symmetry down to SU(L)y. 
Then we have two non-zero-mass states: let us discuss them. In Section 14 (after having 
derived some useful expressions for the chiral condensate) we will see that, in the limit of 
large number of colours N c , we have the following iV c -dependences: 

F w = 0(N^ 2 ); F x = 0(N 1 J 2 ); A = 0(1); c = 0(N B ). (11.12) 



The two eigenstates 0i and 02 corresponding to A^_i and Xl are, in the N c — > oo limit: 

h = , 1 (VLF x Sx + F w S w ), 
'F 2 + LF 2 



i 1 (-FA + VlFxS^), (11.13) 



where SV is the usual "quark-anti-quark" SU(L) -singlet field (in terms of the fields 7Tj it 
is given by: S w = ^Tr$ = A^Y^Ui^i}- So Al-i is the squared mass of the field 
0i and Xl is the squared mass m| 2 of the field 02- At the leading order in the N c — ► oo 
limit, they are given by: 



2LA 



1 _ " fz + lfV 



7T- 1 X 



Let us observe that is of order C(-^), while m| 2 is of order 0(1). As a check, we 
immediately see that, if we put Fx = in the above formulae (i.e., if we neglect the new 
chiral condensate), then 0i = SV and reduces to ^=^, which is the "usual" r)' mass in 
the chiral limit || [TIJ. Moreover m| 2 goes to infinity when Fx goes to zero, so that the 
field 02 = Sx is forced to be zero. Yet, in the general case Fx ^ 0, the two states which 
diagonalize the squared mass matrix A are linear combinations of the "quark-anti-quark" 
singlet field S n and of the "exotic" 2L-fermion field Sx- One of them (0 X ) has a "light" 
mass, in the sense of the N c — > oo limit, being = O(jj-); this mass is intimately 
related to the anomaly and they both vanish in the N c —>■ oo limit (see the first Eq. 
( |11.14| )). On the contrary the field 02 has a sort of heavy "hadronic" mass of order 0(1) 
in the large- N c limit (see Ref. [[JIJ for a detailed discussion about hadrons in the 1/N C 
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expansion). Both the 0i and the 02 have the same quantum numbers (spin, parity and so 
on), but they have a different quark content: one is mostly SV ~ iJ2f =1 {quint — QirQil), 
the other is mostly S x ~ i[det(q sL q tR ) - det(q sR q tL )}. 

What happens when approaching the chiral transition temperature T c ^? We know 

2A 



that F^iT) — > when T — > T c h. From the first Eq. ( |11.14j ) we see that m^T, 



'-ch 



which is just the squared mass m| x of the field Sx for T > T ch : and in fact, from the 
first Eq. (|11.13|) , 0i(T c /j) = Sx- We have continuity in the mass spectrum of the theory 



through the chiral phase transition at T = T c ^. If we assume that c(T) goes to zero at 
least as F%(T) when T — > T c h, then also m^ 2 has a finite limit when T — > T^: m^ 2 — > ^0. 

Another relevant comment has to be made about the field 0i and its squared-mass 
formula, given by the first Eq. ( |11.14| ). It turns out that 4>i is J us t the meson state, with 



a mass squared of order 1/N C , whose contribution to the full topological susceptibility x 
exactly cancels out (in the chiral limit of massless quarks!) the pure gauge part A of x, 
so making x — 0: this is the so-called Witten's mechanism, which is explained in Section 
4. Let us see how this picture comes out in our theory. Let us determine, first, the U(l) 
axial current, starting from our Lagrangian (|9.15|) . This is easily done, remembering how 
the fields U and X transform under a U(l) chiral transformation (see Eqs. ( |9.7|) and 
( |9.10| )). We thus find the following expression for the U(l) axial current J^'- 



42 = i[rr(rfdpU - Ud^) + Lpft^x _ XdpX% (11.15) 



After having inserted here the expressions ( |11.3| ) and ( |11.5|) in place of U and X respec- 



tively, the current takes the following form: 



= -V2Ld,(F T S w + VLFxSx). (11.16) 

In the case in which Sx = (i.e., if we only have the usual chiral condensate (qq)), Eq. 
11.16|) reduces to 42 = "V^LFs-c^SV, where F$ = F w is the well-known expression 



| |10| for the singlet (S n ) decay constant, at the leading order in the 1/N C expansion. 



Instead, in the general case we are considering, the first Eq. ( |11.13| ) allows us to write 
Eq. ( pTTTBD as: 

42 = -V^LF^l (11.17) 



t(L) 



(where F^ is defined below, in Eq. ( |11.18 )). This means that the axial current J$J only 



couples to the ll lighf field <pi, and not to the "heavy" field 02 ( u light' and "heavy" are in 
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the sense of iV c -order) . The relative coupling between and 0i, i.e., the singlet (0i) 
decay constant defined as (0| ^5^(O)I0i(p)) = ^V^Lp^F^, is now given by: 



F* = y/F* + LF%. (11.18) 

Let us now recall the Witten's argument, which is discussed in Section 4, and write the 
two-point function (at four-momentum q) of the topological charge density Q(x) as a 
sum over one-hadron poles, i.e., one-hadron intermediate states: 



X (q) = -i fd 4 xe^(TQ(x)Q(0)) = A (q) + £ (n.ig) 

mesons 1 " u v 



where A (q) is the pure Yang-Mills contributions from the glueball intermediate states 
and it is the leading-order term in 1/N C . In the chiral limit in which we have L massless 
quarks, the full topological susceptibility x = x(Q = 0) must vanish, for the reasons 
explained in Sections 4 and 5: so there has to be a meson state, with squared mass 
m 2 n = 0{1/N C ) (since A o (0) = 0(1), while |(0|Q|n)| 2 = 0{1/N C )), which exactly cancels 
out A (0). This is the meson we usually call rf\ the other meson states in ( |11.19| ) have 
squared masses of order 0(1), so that their contributions to the summation in ( |11.19|) are 
suppressed by a factor of 1/N C . Therefore we obtain that: 

'<°'W'I 2 =J 4, (11.20) 

where A = Aq(0) is the pure Yang-Mills topological susceptibility in the large-iV c limit. 
It is well-known that, in the chiral limit of massless quarks, the topological charge density 
Q(x) is directly connected with the four-divergence of the axial current J^, via the 
anomaly equation: 

d*jlft{x)=2LQ{x), (11.21) 
so that (0\Q\r]') = -J=m^,F n ', which can be substituted into Eq. (|11.20|) to give: 



m 2 F 2 

A=^t. (11.22) 

Eq. ( |11.22| ) relates the mass rn v i of the rf state, its decay constant F^ and the pure-gauge 
topological susceptibility A. It is easy to see, now, that the <fti state, whose mass and 
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11.14|) and (|11.18|) , really satisfies Eq. 



decay constant are given respectively by Eqs. 
(|Tl~22l) : it is just the rf state! 

One could have also followed an inverse reasoning, starting from Witten's assumption 
that there has to be a meson state, called r/', with squared mass nil, = 0(1/ N c ), which 



couples to in such a way to exactly cancel A = Aq(0) in Eq. ( |11.19| ). Looking at 
Eq. ( 11.16 ) and comparing it with the first Eq. ( 11.13 ), it seems natural to conclude that 
(doing also a suitable state-normalization) the state rf is nothing but 0i, with a decay 
constant given by Eq. ( |11.18|) . One thus derives the 0i mass from Eq. ( |11.22|) , which 
obviously furnishes the usual value ( |11.14| ): = = F 2+l F z ■ Therefore our theory 
seems to be self-consistent. 



12. The real- world case at T < T, 



ch 



Let us apply, now, the results of the previous Section to the 11 real-world 1 case, in which 
there are L = 3 light flavours, named u, d and s, with masses m u , m^, m s small compared 



to the QCD mass-scale Aqcd (or better: m u , <C m s ^ Aqcd) [ill. It is useful to 
represent $ as $ = J2a=i ^aTa + T a (a = 1, . . . , 8) being the generators of the algebra 
of SU(3) in the fundamental representation, with normalization: 

Tr(r a ) = ; Tr (r a r b ) = S ab . (12.1) 

(In other words r a = XJ a/2, where A a are the Gell-Mann matrices.) The fields 7r a describe 
the mesons of the octet, while SV is an S'f/(3)-singlet field. In terms of $ and Sx, the 
quadratic part of the Lagrangian ( |9.15|) reads: 



+ 

2 



If we introduce in place of $ its expression in terms of the fields n a and S n , we immediately 
find that the states 7Ti, 7r 2 , 7r 4 , 7r 5 , 7r 6 , tt 7 are already diagonal, with masses: 

< 2 = ml ± =B(m u + m d ), 
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ni 



ni 



""4,5 
2 



= m 



= m 



K° 



B(m u + m s ), 
R0 = B(m d + m s 



12.3) 



where we have put B 



H?-. On the contrary the states 7r 3 , tt 8 , S n , Sx mix together. 

(m -m ) 

However, if we neglect terms of order 0( a^^ /) as ^ * s usua Uy done in chiral pertur- 
bation theory (they are effects of SU(2) isospin breaking, which are experimentally very 
small), we find that also the state 713 becomes "diagonal", with mass: 



ni 



""3 



B{m u + rrid) = 2Bm, 



12.4) 



m u +m d 
2 



where we have put: m 

The interesting results are obviously obtained when considering the system of fields 
7r 8 , S n , Sx, having the following squared mass matrix: 



A 



fh+2m s 

2B^{m 



V 







2B^(m-m s 

6(Aw?+c) , 

— p2 r Jo 

2Au) 1 (1-uji)-2c 

Kfx 





2Acji(1-oi 1 )-2c 

F*F X 
2A(l-wi) 2 +2c 

w 1 

x 



\ 



12.5) 



where f = ^B{m u + m d + m s ) = |£>(2m + m s ). The eigenvalues of A may be quite easily 
obtained at the first order in mi and 1/N C (so that, for example, terms of order 0(rrii/N c ) 
are considered of second order). For this purpose we need to know the iV c -dependences 
of the various quantities appearing in the matrix A. They will be derived in Section 14, 
with the result: 



F n = 0(Ny 2 ); F x = 0(N^y } A = 0(1); c = 0(N c ). 



1/2 > 



12.6) 



The eigenvalues of A are found to be: 

m + 2m s 



2B- 



6A 



+ 



m 



'IX 



2^ + 3Fl) +0 1 



-B(2m + m s ), 



3F X 



Ft Fx 

7T A 



Nj F? + 3F X 3 



-B(2m + m s 



12.7) 



where Oi(l/N c ) is a quantity of order l/N c , which also depends on the parameter uj\ 
appearing in the Lagrangian (|12.2 ). The physical interpretation of these three states is 
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clear. The 77-state is the eighth pseudo-Goldstone meson of the octet: its mass vanishes 
with the (light) quark masses. On the contrary the r/-state and the r/x-state have masses 
which do not vanish when the (light) quark masses are put equal to zero. Yet, while the 
77'— state has a "non-chiral" topological mass p?^pi , which is of order 0(1/N C ) in the 
1/N C expansion, the 77x-state has a sort of heavy "hadronic" mass of order 0(1) in the 
large- A/" c limit. Both the rj' and the rjx have the same quantum numbers (spin, parity and 
so on), but they have a different quark content: one is mostly ~ iYn=i (QiLQiR — QirQil), 
the other is mostly ~ i[det(q S LQtR) — det(q s Rq t L))- From Eqs. (|12.3|) , (|12.4T) and ( |12.7|) we 
find that the squared masses of the octet mesons satisfy the Gell-Mann-Okubo (GMO) 
formula, which is a standard result of chiral perturbation theory at the lowest order: 

3m 2 v + ml = 4:m 2 K , (12.8) 

considering: m u ~ ~ m. In fact it is natural to suppose that the introduction of 
a new chiral order parameter, which breaks only the U(l) axial symmetry, should not 
modify the mass relations (such as the ( |12.8| )) of the SU(3) meson octet: these relations 
only derive from the breaking of 577(3) <g> SU(3) down to SU(3)v- 

The new interesting result is obtained considering also the second equation of ( |12.7| ). 
From this we immediately derive the following expression 




QA 

m ;,+m;-2m^ = - (12.9) 

7T 

Eq. (|12.9|) is an obvious and natural generalization of the Witten-Veneziano formula for 



the T]', derived in Refs. 0, |10|- Now the presence of the new U(l) chiral order parameter 
induces a correction of order F x /F 2 , where Fx is essentially the magnitude of this chiral 
condensate. It is a beautiful and interesting thing that Eq. ( |12.9j ) does not contain the 
other (unknown) parameters of our model (such as ci, Ui, . . .). One may ask if this is 
true only for our simplified model. In the next Section we will see that the result ( |12.9| ) is 



really "model-independent" , being still valid when one adds to our "simple" model other 
possible additional terms. 
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13. The effects of additional terms in the Lagrangian 



As we have said in Section 9, the form ( |9.13| )-( |9.14 ) for the chiral effective Lagrangian, 
at the leading order in rrti and 1/N C , is surely the most simple and natural one, but not 
the most general one. Now we want to see if our result (|12.9|) is "model-independent", 
and we will find that this is just the case. 

To prove this, we first consider the effects produced by an additional term, linear in 
the quark mass matrix M, of the form: 



5C {2) 



Pi 



Tr(MU) det(f/)X t + Tr(M t f7 t ) det(t/ f )X 



+ 



+P2 



Tr(MU) det(W)X + Ti(M^U r ) det(f/)X t 



:i3.r 



This new term may also be written in the following form: 



s 



Tr(MU) + Tr(MW)] ■ \det(U)X* + det(f/ f )X 



+D \Tr(MU) - Tt{M^)] ■ fdet(C/)X t - det(f/ t )X 



'13.2) 



where the new parameters S and D are linked to p\ and pi by the following relation: 



S 



Pi + P2 



D 



Pi ~ P2 



(13.3) 



Proceeding as in the previous Section, we find that the squared masses of the pseudoscalar 
mesons are given by (expanding to the first order in and 1/N C ): 



m 
m 
in 



2 

7T1,2,3 
2 

•""4,5 
2 

""6,7 

ml 



m r 



in 



'IX 



ml = B f (m u + m d ), 
m 2 K± = B'(m u + m s ), 
m 2 R0 = B>(m d + m s 



IB 



,m + 2m s 



6A 



+ 



F2 + 3F1 + 3 



F 2 2 

-B'(2m + m s 



x 



2c'(F 2 + 3F 2 ) 1 
FjF 2 + 0l{ N~, ) + 



3F 2 X 



F2 + 3F 2 3 



2V2F 2 F x D(2m + 



B' (2m + m s 



13.4) 
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where the new parameters B' and d are given by: 



B' 



B + ^F X S, 
c+^-F*F x S(2m + m s 



13.5) 



Observing Eqs. ( |13.4| ), we immediately see that both the GMO relation ( 12.8 ) and our 



generalized Witten-Veneziano formula ( |12.9| ) are still valid. The effect of the introduction 
of the additional term ( |13.2[) to the effective Lagrangian (|9.13 )- (|9.14 ) on the masses of 
the octet mesons and of the rj' is therefore only that of rescaling the parameter B — > 
B' = B + ^F^FxS (note that only the parameter S appears here): the relations (|12.8|) 
and ( |12.9| ) do not depend (in form) on the value of B, so that they remain valid. In a 
certain sense, the "model-dependence" is absorbed entirely by the squared mass rn^ x of 
the 7/x-state, which is "large" in the sense of the 1/N C expansion, being of order 0(1). 

By a similar analysis, it is not hard to verify that also the introduction of other (more 
complicated) additional terms, linear in M, has the same effect: it simply rescales 5* and 
D in Eqs. ( p^4] ), and the relations flETgp and ([T2T9| ) are still valid. 

One may also try to generalize the last term in the potential (|9.14|) by adding, for 
example, a term of the form: 



5V 



C-2 



8V2 



det(f/)X t + det(f/ t )X 



13.6) 



It is immediate to see that its only effect on the quadratic part £2 of the Lagrangian 
written in Eq. ( p,2.2| ) is to re-scale the parameter c: 



V2\V2j \V2, 



(13.7) 



One can easily convince oneself that a general modification of the above-mentioned term 
in the potential V(U, U\ X, X*) (Eq. (|13.6|) is an example of such a modification) implies 
only a re-scaling of the parameter c in the quadratic Lagrangian ( |12.2| ). Once again all 
the "model-dependence" is absorbed by the squared mass m^ x of the 7/x-state. 



The GMO formula ( |12.8| ) and the generalized Witten-Veneziano formula ( |12.9| ) acquire 
the aspect of theorems: they are not sensitive to the specificity of the model. And while 
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the GMO formula does not depend on the new U(l) chiral order parameter (as it must 
be!), we find that the Witten-Veneziano formula (which links together the rj mass, the 
77 and K masses and the pure-gauge topological susceptibility A) is modified through 
a term which only contains Fx, i.e., essentially the magnitude of this U(l) chiral order 
parameter. Experimentally we have that: 

m K± = 493.646 ± 0.009 MeV, 

m KOjRO = 497.671 ± 0.0031 MeV, 

m v = 547.45 ±0.19 MeV, 

m v , {958) = 95 7.75 ±0.14 MeV. (13.8) 

From lattice simulations for the SU(3) Yang-Mills theory (QCD with no quarks) one 
obtains the following value for the topological susceptibility A (see, e.g., Refs. ||| 
and pi): 

A ~ (179 ± 4 MeV) 4 . (13.9) 



Then Eq. ( |12.9| ) implies that, within the present errors, a value for Fx different from 



zero is not excluded, up to a superior limit of about 10 -j- 20 MeV (we remind, for a 
comparison, that the pion decay constant is about F n ~ 93 MeV). Of course, one should 
also remark that Eq. ( |12.9| ) has been derived in the large-iV c limit, and possible next-to- 
leading corrections in the 1/N C expansion could be non-negligible at the physical value 
N c = 3 (see, for example, Ref. ]50|]), so becoming comparable with (or even larger than) 
the (presumably!) small correction ~ F x /F%, induced by the new U(l) chiral condensate. 



14. The topological susceptibility and the chiral condensate 

After having discussed in the previous Sections the mass spectrum of the theory above 
and below the chiral phase transition at T = T c h, we can now see which are the predic- 
tions of our model about the topological susceptibility and the chiral condensate. Let 
us begin with the topological susceptibility, defined as (TQ(x)Q(0)) FT -(k = 0), where 
"F.T." stands for the Fourier transformed. If we want to derive this two-point function 
of Q(x), we need to consider the Lagrangian in the form ( |9.13| ), where the field variable 
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Q(x) has not yet been integrated. So doing, we obtain the following expression for the 
two-point function of Q(x), in the region of temperatures T ch < T < Tuny. 



(TQ(x)Q(0)) RT ik) =*(A-\k)) u 



i4.r 



where A 1 (k) is the inverse of the following squared mass matrix for the ensemble of fields 
(Q(x), Sx, bn, b\2, ■ ■ •) (see Section 10): 



/ i 



A(*) 



V2 
Fx 




V2 
Fx 



0(m 





0{m L - 1 ) 
k 2 - \\lBl 

2 7T 7T 



V 



\ 



/ 



14.2) 



where m\ is the eigenvalue proportional to det(M) = rife=i m fc, written explicitly in Eq. 
( 10.3|) . Doing explicitly the calculations, we obtain that: 



(TQ(x)Q(0)) RT ik) = i(A- l (*)) u = iA k ~ 



k 2 - m 2 Sx ' 



14.3) 



where, as usual: tn 2 Sx 



Therefore the topological susceptibility in this region 



T ch <T < Tu(i), is given by: 



(TQ(x)Q(0)) RT ik = 0) = L4l^_ = iA 



m 



Sx 



Fx 2 1 



sup( mi )-0 2 2 \yf2\lBl, 



det(M). 



14.4) 



So we obtain that the topological susceptibility of the full theory (not of the pure YM 
theory!) in the region T c h < T < TJ/fi) is proportional to the product of the light quark 
masses rrii (i — 1, 2, . . . , L). This is the same behaviour that can be derived from a simple 
instantonic model, as shown in the Appendix. We also find that the propagator of the 
field Sx is given by: 



(TS x (x)Sx(0)) F - T ik) = i (A-\k)) 



22 



k 2 — rn 2 Sx ' 



(14.5) 
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which is consistent with what was found in Section 10: i.e., the field Sx has a squared 



mass m% = mi + -§#. 

bx F x 

What happens in the region T < TdJ Here we have to consider the following quadratic 
Lagrangian, including the field Q(x) (see Sections 9 and 11): 

A = ^E + ^ E d^d^ + ±d,S x d»S x + 



-oLW "oL o 77 77 + 



Z i=l Z 



1 fV2J^„ V2„ \ 2 1 



— -c 



2 F,^ F x A 2A 



i=l 

-oo x ^Q{x) E vfi - (1 - ^)^g(x)5 x . (14.6) 

Proceeding as in the previous case, we derive the following expression (in the special case 
L = 2) for the topological susceptibility in the region T < T ch , below the chiral phase 
transition: 

2 2 

(TQ(x)Q(0)) RT ik = 0) = iA—^ ™- , (14.7) 

c ^1^2 + VP1 "T A*2^ 

which has the usual dependence on the light quark masses in this region: it goes as 
oc m\ for m\ — > (Using Eq. ( |14.9| ), that we will derive below, and remembering 
that /i z 2 = fpmj, we derive that (TQ(x)Q(0)) F - T -(k = 0) = -imi(qiqi), for T < T c/l and 
— > 0, as requested by the QCD Ward Identities at the leading order in the light quark 

0)- " 



masses 



Now we address the question of the chiral condensate. It is well known that the deriva- 
tive of the QCD Hamiltonian with respect to m« is the operator q^ (being SCqqd = 
— J2i=i m iQiQi) ■ The corresponding derivative of the vacuum energy therefore represents 
the vacuum expectation value of q\qi. So, in our case, it must be: 

(<hQi) = -^ i (V(U,rt,X,Xt)), (148) 

where V(U,lP,X,X*) is the potential term written in Eq. (|9.14 ). It is immediate to 
calculate the expectation value of the potential V in the two regions T < T c h and T c h < 
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T < Tuny In the region T < T c h we find that (always at the leading order in mj, 1/A|, 
1/A 2 x): 

fe>T<T e „ ^ -^-Bm^Tr- (14.9) 

It is of order one in the light quark masses, indicating that the SU(L) <g> SU(L) chiral 
symmetry is really broken. Of course the interesting result is obtained in the region 
T ch < T < Tum, above the SU(L) ® SU(L) chiral transition. Here we find that: 

(w*)T ch <T<T u{1) ^ -t?|^ - \c x F x (-^-] (II m k ). (14.10) 



The interpretation of the two terms in the right-hand side of Eq. ( 14.10| ) is rather simple 



in a diagrammatic language. The first term, linear in the mass mj, corresponds to a 
diagram with a mass insertion, as shown in Fig. 5. 




Figure 5: The first term in the RHS of Eq. Q14.10Q . 

Instead the second term clearly corresponds to an insertion of the new 2L-fermion effective 
vertex associated with X, as shown in Fig. 6. 

In conclusion we have obtained that, in the region T c h < T < Ty-m, the chiral condensate 
(qiqi) has the following behaviour with respect to the light quark masses: 

(q i q i ) = 1 {m i )+0 2 {l[m k ). (14.11) 

This relation can be derived also in a simple instantonic model, as shown in the Appendix. 
Comparing the explicit expression for 02{J\k^i rn k)i written in the right-hand side of Eq. 
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Figure 6: The second term in the RHS of Eq. (|14.10|) . 



( |14.10| ), with the expression for the full topological susceptibility given by Eq. ( |14.4| ), we 



find the following interesting relation: 

(TQ(x)Q(0)) KT ik = 0) = -imi0 2 ([[ m k ), (14.12) 

which is, of course, independent of the flavour index "Z" . 

We want to close this Section with some comments about the ^-dependence of some 
relevant quantities that we have considered in our analysis. In Section 8 we have shown 
that the picture we have discussed so far is really consistent with the saturation of the QCD 
Ward Identities above the SU(L) ® SU(L) chiral phase transition [^j. The saturation is 



obtained in the framework of the 1/N C expansion, where N c is the number of colours and 
we derive the N c dependence of a certain number of relevant quantities. In particular the 
squared mass m 2 Sx of the pseudoscalar singlet boson in the region T ch < T < TJy^), and 
in the chiral limit of zero quark masses, depends on N c as Q(l/N c ). From Eq. ( |10.5| ) 
we have that m| = J4- for zero quark masses; considering that A = 0(1) with respect 
to N c (see Refs. | 0), we immediately derive that F x = 0{N l J 2 ). In Eq. (gAg ) we 



have found that the chiral condensate in the region T < T c h is given by: (qiqi)T<T ch — 
— ■^B m F n . It is well known (see Refs. U and [|^, j3^| ) that both F 2 and (qiqi) are of 



order 0(N C ) (and this remains true also considering the new condensate which breaks 
the U(1)a symmetry). This fact implies that B m = 0(N^ 2 ). In Section 8 we have also 
found that the two contributions Ci(mj) and 2 (J\k^i m k) for the chiral condensate (qiqi) 
in the region T c h < T < Tjjn\ (see Fig. 4 in Section 14) are both of order 0(N C ). This 
means, considering the explicit expressions for Oi(rrii) and C^dlfc^i^fc) written in Eq. 
Piq ), that \\Bl = 0(1) with respect to iV c and that c x = 0{N^/ 2 ). Therefore, 
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from Eq. ( [LI .8 ), c = 0(N C ). As a consequence of this, we have derived in Section 11 the 
iV c -dependence of the non- zero-mass states <f>\ and 02 (whose explicit expressions in terms 
of the fields S n and Sx are given in Eq. ( |11.13| )) in the region T < T ch . We have found 
that the state <fii has a "light" mass m^, as JV C oo, being = p^+lf' 2 = ®(jt)- 
Moreover is purely originated by the anomaly (it vanishes if A is put equal to zero!). 
On the contrary the state (f>2 has a sort of heavy "hadronic" mass m^ 2 of order 0(1) with 
respect to N c and not dependent on the anomaly. 



15. Conclusions 

In this paper (after a pedagogical review of the U(l) problem at T = 0, which has been 
presented in Sections 2-5), we have tried to gain a physical insight into the breaking 
mechanism of the U(l) axial symmetry, through a study of the behaviour of the theory 
at finite temperature. As discussed in the Introduction, the topological susceptibility A 
of the pure Yang-Mills theory is a fundamental quantity for studying the U(l) chiral 
symmetry, both at zero and non-zero temperature. From some previous works of Witten 



and Veneziano [pl| |16j, it is known that a value for A different from zero implies, at 
large number N c of colours, the breaking of the U(l) axial symmetry, since it implies 
the existence of a pseudo-Goldstone particle with the same quantum numbers of the rf . 
In the Witten- Veneziano model the mass ms v of the SU (L)-singlet field for T < T c h 
acquires a "topological" contribution from the topological susceptibility in such a way 
that, in the chiral limit of zero quark masses, m| = ^M- (F n is the usual pion decay 
constant: F v ~ 93 MeV at T = 0). Approaching the chiral phase transition temperature 
T c h, F n vanishes and so we could expect that also A goes to zero approaching T c h, in 
order to avoid a singularity in the singlet field mass for T —* T c h- Yet, recent lattice 
results [|nj show that the YM topological susceptibility A(T) is approximately constant 
up to the critical temperature T ch , it has a sharp decrease above the transition, but it 
remains different from zero up to ~ 1.2 T c h- Also present lattice data for the so-called 
"chiral susceptibilities" (discussed in the Introduction) seem to indicate that the U(l) 
axial symmetry is still broken above T c h, up to ~ 1.2 T c ^. 

In the following, we briefly summarize the main points that we have discussed and the 
results that we have obtained. 
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• One expects that, above a certain critical temperature Tu^, also the U(l) axial 
symmetry will be (approximately) restored. We have tried to see if this transition 
has (or has not) anything to do with the usual SU (L) ® SU (L) chiral transition: 
various possible scenarios have been discussed in Section 6. 

• We have proposed a scenario (supported by the above-mentioned lattice results) 
in which the U(l) axial symmetry is still broken above the chiral transition and 
the pure YM topological susceptibility vanishes at a temperature T x between T ch 
and Tu(i). A new order parameter has been introduced in Section 7 for the U(l) 
axial symmetry. In Section 8 we have shown that this picture is consistent with 
the saturation of the QCD Ward Identities above the SU (L) <g> SU (L) chiral phase 
transition, in the framework of the 1/N C expansion, taking care also of the new 
2L-fermion condensate which breaks the U(1)a symmetry. 

• In Sections 9-14 we have shown that this scenario can be consistently reproduced 
using an effective Lagrangian model. We have analysed the effects that one should 
observe on the mass spectrum of the theory, both below and above T ch . In particu- 
lar, the Witten-Veneziano formula for the mass of the rf is modified by the presence 
of the new 2L-fermion condensate. In this scenario, the r{ survives across the chi- 
ral transition at T c h in the form of an "exotic" 2L-fermion state ~ i[det(q S LqtR) 
— det(q S Rq t L)]. This particle is nothing but the would-be Goldstone boson coming 
from the breaking of the U(l) axial symmetry. For T > T c h, it acquires a "topo- 
logical" squared mass of the form 2A/Fx, where F x is essentially the magnitude of 
the new order parameter for the U(l) chiral symmetry alone. 

This scenario could perhaps be verified in the near future by heavy-ions experiments, by 
analysing the spectrum in the singlet sector. If the "exotic" light state would be copiously 
produced at high T, its subsequent decay into L light mesons (as the temperature T of 
the plasma formed in the collision decreases), could be seen in terms of a production of L 
quark-antiquark couples (i — 1, . . . , L) and their subsequent linear superpositions to 
form L physical mesons. This is quite similar to high energy instanton effects. Some tests 
and verifications of this picture could also be provided by Monte Carlo simulations on the 
lattice. We hope that further progress along this line will be done in the near future. 
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Appendix: the quark mass dependence for T > T c h 



In this Appendix we derive a (possible) quark mass dependence of the topological suscep- 
tibility x °f the full theory (QCD with quarks) and of the chiral condensate {qiqi) in the 



chirally symmetric phase (i.e., for T > T c h). We follow a procedure outlined in Ref. [51 
and we work in the Euclidean theory. It is well known that: 

lim = -vrp(O), (A.l) 

sup(m fe )^0 

where p(A) is simply the gluonic expectation value of p(X;A), which is the density of 
eigenvalues of D[A] = (-D M = + igA^). In the phase where chiral symmetry is 

restored we have that: 

lim (q iqi ) T>Tch = 0, (A.2) 

sup(m ft )^0 

so that p(0) = 0: i.e., the "average" density of eigenvalues of D[A] near A = is zero 
in the chirally symmetric phase. Recall that the partition function becomes (Sq is the 
gluonic action): 

Nf 

Z = f [dA}e~ SG ^ J] det(D[A] + M k ), (A.3) 
J k=i 

if we integrate over the Nf quark fields (L of which are taken to be "light" quarks of mass 
m, Kqcd-i Aqcd being the QCD mass-scale; the chiral limit is done for this L flavours: 
sup(mj) — > 0). If we are in the chirally symmetric phase, then any discrete zero mode 
due to a non-zero topological charge Q of the gauge field A will be well separated from 
the continuous distribution of p(\), since this latter distribution goes to zero at A = 
by virtue of (|A.1|) and ( |A.2| ). The Atiyah-Singer index theorem relates the number of 



left-handed n_ and right-handed n + zero modes of the fermionic operator D[A] to the 
topological charge Q of the background gauge field: 

n+ -n_ = Q. (A.4) 

So, if we define ti[-A]q = n + + n_ — \Q\ [n + + n_ being the total number of discrete zero 
modes), we find that: 

oo /N f \ Q+ n l A h 

Z =H! [dA] Q e- s °W« [YImA f(M, [A]q) = 

Q=0 J \i=l J 
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oo /'Nf \^ /Nf \ "M« 



E II M A J [dA] Q e- SG[A] ^ (UmA f(M, [A]q 



=0 \i=l 
x /'Nf 



where: 



e n^i «qm, (a.5) 

Q=0 \i=l 



/N f \ n i A h 

a Q (M) = J [dA] Q e- s °W* I^mA f(M, [A] Q ). (A.6) 



The f(M, [A]q) picks up the contribution of the eigenvalues A different from zero: it is of 
order 0(1) in the chiral limit sup(mj) — > 0. Moreover, it is always possible (at least for 
Q = 0, Q — 1) to find a configuration of the gauge field A in which = 0. All of 

this implies that aq(M) is regular and, in general, non-zero for sup(mj) — > 0. If we are 
in a fixed physical four-volume V4, then we clearly have: 

(Nf \ Q 

T,Q 2 [ll M A «q( m ) 

1 , 2* Q=0 \i=l 



V 4 00 (Nf \ Q sup(m fc )-»0 

^e (ri M ®q( m ) 

=0 \i=l / 




Mi ^(0 0^,...,^,) = o(det(M))) (A7) 



fc=L+l 



7 4 ao(0,...,0 J M L+1 ,...,M A r / ; 



where M is the usual mass matrix of the L light quarks (m^ = M, for z = 1, 2, . . . , L). That 
is, x = C(lli=i m i) m ^ ne chirally symmetric phase. The uncertainties of this method 
lie in the finite-volume approximation and in the WKB approximation: probably this 
is valid here owing to the fact that p(0) = 0. Viceversa, in the region where the chiral 
symmetry is broken we have that p(0) 7^ 0, so that any discrete zero mode is embedded in 
a "sea" of continuous zero modes. That's why it is not immediate to extend the method 
used before also to the case T < T ch . In Ref. the ( |A.7| ) is derived in a different way, 



without the WKB approximation, yet always in a finite volume. With the same technique 
we can also calculate the chiral condensate (ftft). It is clear that: 

N f 

-S a [A)-J d±x Qki^D^ + M k )q k 
,- v _ f[dA] [dgpgfae k=i 

\QiQi) ry 
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J d_ 

V A ZdMi 
1 dZ 
ViZdM; 



[dA] [dq] [dq]e 



k=l 



S G [A]- J d^x QkirfD^ + M k )q k 



(A.8) 



where V4 is always the physical four-volume. Inserting here the (|A.5 ), we find that: 

(QiQi) z 



I r) 00 / Nf 

E II^I a Q (M) 



V 4 Z 



V 4 Z dMi 



1 00 

E 

Q=0 



Q=0 \k=l 

k — 1 J ^ \ k^^i 



(A.9) 



In the chiral limit sup (m;) — > 0, the (|A.9| ) becomes: 



sup(m fc )^0 

<9a 
dMi + 



Y[M k )a 1 (0,...,0,M L+u ...,M J 



X 



V 4 ao(0,...,0,Afi + i,...,Mjv > ; 



(A.10) 



Considering that, if the operator -D[A] has the eigenvalue A, it also has the opposite 
eigenvalue —A, one immediately derives from the ( |A.6| ) that = 0(rrii). And so we 
have that: 

(m) = 1 (m i ) + 2 (l[m k ). (A.ll) 

k^i 



The results ( |A.7| ) and ( A. 1 1| ) are the same that we have found from our model in the region 
T c h < T < Tum (see Section 14). Moreover, from the explicit expression of C^dlfc^i^fc) 
given in Eq. ( |A.10| ), we derive, after a comparison with Eq. ( |A.7| ), the following relation: 



-(Q 2 ) = -m^dlrrik) 



(A.12) 



(which is, of course, independent of the flavour index "i"). This is just the Euclidean 
version of Eq. (|14.12|) , derived from our effective Lagrangian model. 
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